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Abstract 

> 

00 ' We study zero modes of two-dimensional Pauli operators with Aharonov-Bohm 

■ fluxes in the case when the solenoids are arranged in periodic structures like 

, chains or lattices. We also consider perturbations to such periodic systems which 

may be infinite and irregular but they are always supposed to be sufficiently 

' 

o 

1. Introduction 

^ ' The appearance of zero modes (wave functions at zero energy which are ground states 
for a positive quantum Hamiltonian) belongs to the most interesting phenomena in 
systems with topologically non-trivial configuration spaces; see the discussion and an 
^ . extensive bibliography in Zero modes of the Dirac and Pauli operators are of 
great importance in many places in quantum field theory and mathematical physics 
Oinilll. They are the ingredients for the computation of the index of these operators 
and play a key role in understanding anomalies. One of the best known examples for 
such operators is the Pauli Hamiltonian of a two-dimensional charged particle moving 
in a magnetic field perpendicular to a plane and penetrating the plane in a bounded 
domain. In this case the field defines a vector bundle with a non-trivial connection 
and zero modes appear at sufficiently high strength of the field jSj. More precisely, it 
is easy to prove that the dimension d of the space of zero modes is d = [|$|J where 
$ is the total flux of the magnetic field measured in magnetic flux quanta, and for a 
real x, x > 0, [xj denotes the lower integer part of x ([Oj =0, \n\ = n — 1 for n > 1 
integer, and otherwise [xj = [x], the integer part of x). It is worthy to note that in the 
three-dimensional case the appearance and the degeneracy of zero modes is a more 
subtle fact (see e.g. [El 13 El El CHI IH] and the discussion therein). 



1 



2 



In the current paper we restrict our consideration to two-dimensional systems only. 
More precisely, we consider Pauli operators which are Hamiltonians of an electron 
confined to a plane and subjected to a perpendicular time-independent magnetic field 
which is the sum of a uniform field and an additional field contributed by a (finite 
or infinite) array of singular fiux tubes or, in other words, by an array of solenoids 
of zero width. We focus on zero modes in such systems. In more detail, the aim of 
the paper is to find conditions for appearance of zero modes in systems placed in a 
magnetic field with an infinite array of Aharonov-Bohm vortices. It has been shown 
in ^2] on the physical level of rigor that zero modes occur if Aharonov-Bohm vortices 
are arranged in a periodic plane lattice provided that not all magnetic fiuxes involved 
have integer values. In this paper we present a rigorous proof and show that under the 
same condition imposed on the fiux, the result is true for a chain of Aharonov-Bohm 
solenoids or, more generally, for a uniformly discrete union of such chains. Moreover, 
the zero modes are retained if one adds to such a periodic structure of Aharonov-Bohm 
solenoids a not necessarily regular array of solenoids having sufficiently low density. 
This stability of zero-modes for the Hamiltonian that we call ifmax (its definition is 
discussed in Section |2} shows that their origin differs from that for localized states in 
the so called Aharonov-Bohm cages ITij , the latter are destroyed by arbitrarily 
small period modulations [THj. 

The main results of the paper are obtained with the help of a version of the 
Aharonov-Casher ansatz 5j. This version was proposed by Dubrovin and Novikov in 
fT^ who employed it for an explicit construction of ground states of periodic magnetic 
Schrodinger operators (see Novikov's review paper [17]). In our case, this ansatz re- 
duces the problem of finding zero-modes to some estimates for entire functions. The 
mechanism of appearance of zero modes in the considered cases is close to that for 
a two-dimensional system in a uniform magnetic field in the presence of an infinite 
array of point scatterers [HHl CHI CHI EDI EH El ■ 

An interesting physical consequence of our result is the occurrence of oscilla- 
tions of the type "localization-delocalization" in periodic systems of Aharonov-Bohm 
solenoids placed in a varying uniform magnetic field (Theorem I8.1(i|l . Another inter- 
esting result described in Theorems 18 . 81 and 18 . 1 61 is related to the problem of absolute 
continuity of the spectrum of the Schrodinger operator with periodic vector poten- 
tial A. This absolute continuity has been proved for a wide class of potentials A 
(23 EH ESI EHl EH EH] • An example of a vector potential A having eigenvalues in the 
spectrum of the corresponding Schrodinger operator was given in [22] but only for 
dimensions higher than 3. Our results give such an example in dimension 2. 

The paper is organized as follows. In Section El we try to point out some aspects re- 
garding the history and the background of the problem. In Section|Hlwe discuss shortly 
the gauge invariance in the case when the magnetic field is a distribution. In Section 0] 
we introduce several basic examples of models with Aharonov-Bohm fiuxes some of 
them are the main subject of this paper and are studied in detail in the sequel. Sectional 
is devoted to a rigorous definition of the Pauli operator with Aharonov-Bohm fiuxes. 
In Section IHl we discuss the elimination of integer- valued Aharonov-Bohm fiuxes. In 
Section 13 we recall the Aharonov-Casher ansatz which makes it possible to construct 
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ground states of the Pauli operator using the theory of analytic functions. The main 
results of the paper are contained in Sections|Hland|ni In Section|Hlwe study zero modes 
of the Pauli operator with an infinite periodic system of Aharonov-Bohm solenoids. In 
Section ini we address the question of perturbations of such periodic structures caused 
by translations and additions of Aharonov-Bohm solenoids. The subsystem formed 
by solenoids affected by the perturbation may be infinite and irregular but we always 
suppose that it is sufficiently scarce. Here we also discuss some examples of irregular 
Aharonov-Bohm systems. For the reader's convenience we have included three ap- 
pendices. In the first appendix we collect some basic definitions and auxiliary results 
concerning lattices. In the second appendix we recall some basic notions and results 
from the theory of analytic functions related to the growth of entire functions. The 
third appendix is devoted to the Weierstrass a-function. 

2. Additional comments on the history and the back- 
ground of the problem 

There are many interesting and important physical problems related to systems in- 
volving Aharonov-Bohm fiuxes. Since the publication of the original paper due to 
Aharonov and Bohm [Hllj the physics of a magnetic fiux in an infinitely thin solenoid 
(called Aharonov-Bohm flux or Aharonov-Bohm vortex) has been investigated both 
from theoretical and experimental points of view • The physical origin of the 

Aharonov-Bohm effect is even a subject of theoretical investigations up to now |33j . 
On the other hand, the motion of a charged particle (an electron, a hole or a composite 
fermion) in a plane perpendicular to a uniform magnetic field has found an important 
application in physics of the quantum Hall effect Ei] • The most striking feature of 
the Hamiltonian of such a system is the Landau quantization of the spectrum which 
consists of highly degenerated equidistant energy levels; this makes quantum Hall phe- 
nomena possible. Moreover, it is of interest to know how the quantum Hall system 
is altered by various defects, in particular, by impurities or by inhomogeneities of 
the magnetic field. Additional Aharonov-Bohm fiuxes appear to be a minimal mod- 
ification of the uniform magnetic field, while general inhomogeneous magnetic fields 
are extremely difficult to handle Similarly, a minimal perturbation of the 

quantum Hall system is given by a point perturbation of the Landau operator (i.e., 
the Schrodinger operator with a uniform magnetic field) [SB]- As shown below, both 
modifications require the operator extension theory for a correct construction of the 
corresponding Hamiltonian |39j . 

The vector potential of a system of Aharonov-Bohm solenoids has a strong singu- 
larity at the points where the plane intersects the solenoids. Therefore the differential 
operator defining the Hamiltonian is not essentially self-adjoint on its natural domain. 
This is true both in the non-relativistic case (for the Schrodinger and Pauli opera- 
tors) and the relativistic one (for the Dirac operator). The boundary conditions for 
Schrodinger operators with an Aharonov-Bohm vortex as well as the corresponding 
self-adjoint extensions (i.e., Hamiltonians describing a spinless non-relativistic quan- 
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turn particle) are considered in many papers, let us mention e.g. ^OJEIESIISI- The 
multi-solenoid case is more difficult because of the rotational symmetry violation. This 
case was treated by means of the Krein resolvent formula in 03], and for an infinite 
chain of solenoids in jU] ; different approaches are presented in |461 ITTl HB] . The prob- 
lem of defining the boundary conditions at the presence of a uniform background field 
has been investigated in jlHl HO] • In the relativistic case, the problem of defining the 
appropriate Dirac operator is discussed e.g. in [HH E21 ISHj, and at the presence of a 
uniform component - in the recent articles |54 [ l55 | IHI) t l57j. In all the mentioned papers, 
the spectral or scattering properties of the derived Hamiltonians are studied as well. 

On sufficiently smooth functions from L^(M^)(8)C^ = L^(]R^; C^) the two-dimensional 
Pauli operator for a charged particle with the spin s and the gyromagnetic ratio g acts 
ClS db formal differential operator |^ 

where e and m* are the charge and the mass of the particle, respectively, A = [A^, Ay) 
is the vector potential of a magnetic field B = Be^, B = dxAy—dyA^, fi is the magnetic 
momentum operator, 

A = gSfiBSz , 

with fiB being the Bohr magneton, fis = — |e|/l/(2m*c), and 

. _ 1 _ 1 / 1 \ 
~ 2 ~ 2 V -1 J 

(we consider the motion of a particle in the plane canonically embedded in the 
space M^). In general, the non- relativistic limit of the Dirac equation leads to the 
value g = 2, and the main part of our work deals with this value of the gyromagnetic 
ratio. In the case of an Aharonov-Bohm solenoid B is proportional to the Dirac 
delta function, S{r), and therefore the operator takes the form of the Schrodinger 
operator perturbed at a point and with a finite coupling constant a standing in front 
of the "5-potential" . On the other hand, it is well known that in the two-dimensional 
case under consideration the expression defines a non-trivial perturbation of the 
operator 

only if a is in some sense "infinitesimal" [S^ (we suppose that appropriate boundary 
conditions defining Hq are chosen). This problem has been analyzed in jSH] in detail 
for an arbitrary positive value of g. To get around it, a solenoid of finite radius R 
is considered with a rotationally symmetric magnetic flux inside the solenoid but 
otherwise having an arbitrary proffie (including the magnetic flux supported on the 
surface of the infinite cyhnder), and the hmit R ^ is discussed. In addition to 
let us also mention papers [^01 ESI 1121 EH EH EH EHl EH EH] • Of course, the same 
approach is useful when a uniform component of the field is present or in the case of 
the Dirac operator (see [HH ] 170 ] ITT] and references therein). 
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In the most important case when the gyromagnetic ratio g equals 2 the Pauh 
operator has remarkable supersymmetry properties which makes it possible to use the 
Aharonov-Casher decomposition j^. As a result, we have a convenient definition of the 
Pauli operator with a singular potential by means of a quadratic form (see Sectional). 
More precisely, in this case we have, as usual, two natural quadratic forms associated 
to the expression (0) - the minimal and the maximal one (with the definition of the 
magnetic Schrodinger operator taken from [Z^). These forms provide us with two 
natural types of self-adjoint operators denoted -ff^ax -ff^in playing the role of 
Pauli operators with Aharonov-Bohm solenoids (the sign ± stands for spin up and 
spin down supersymmetric partners, respectively). There is an important distinction 
between the operators -ffmin and -ffmax- As it follows from the definitions, both operators 
-^min coincide with the Friedrichs extension of the symmetric operator defined by 
expression Q with the vector potential A corresponding to a system of Aharonov- 
Bohm fiuxes. Therefore this extension (denoted simply by ifmin) may be interpreted 
as the Hamiltonian of a "spinless" particle moving in a system of Aharonov-Bohm 
fluxes (this corresponds to physical problems for an electron when the spin-orbit 
coupling can be neglected and spin splitting is taken into account with the help of the 
perturbation theory [SB])- Such a Hamiltonian has been considered e.g., in jmEHI- 
On the other hand, the operators -ffj^ax ^lo not coincide in general which indicates that 
they directly take into account the energy of the spin-orbit interaction and therefore 
they may be regarded as the Pauli operators of the system under consideration. In 
the present article we concentrate mainly on zero modes of -ffmax- Note that boundary 
conditions defining the Hamiltonian -ffmax are given in ^21 HB] (in the case of a single 
solenoid) and in [ZSj (in the two-solenoid case). 

For a finite system of Aharonov-Bohm solenoids, the existence problem of zero- 
energy eigenfunctions was considered in [TH [73 [THl EZl IZHl- In this case the number 
d of linearly independent zero-modes depends on the fractional parts of fiuxes in 
the individual solenoids, {x} = x — [x], rather than only on the total fiux $ in the 
system. This phenomenon is a consequence of the gauge invariance properties for 
the Aharonov-Bohm fiuxes (see e.g. papers [ZOl IHOl IHH IH2I)- In the case when the 
considered magnetic field has a "regular" component in addition to the magnetic field 
of Aharonov-Bohm solenoids the appearance of zero modes has been analyzed in 
|83| IMj . The results of fSl^ are applicable also to the case when an infinite number 
of Aharonov-Bohm solenoids is present in the system but the total magnetic fiux is 
necessarily finite (moreover, after some gauge transformation the total variation of 
the flux must be finite). On the other hand, it is clear that the thermodynamic limit 
of a bounded system with a fixed density of Aharonov-Bohm fiuxes is a system with 
an infinite number of Aharonov-Bohm solenoids and with an infinite total fiux. An 
example for a system of such a kind is the quasi-two-dimensional system with columnar 
defects in a uniform magnetic field directed along the defect axis |HS1 EEl IHI] or the 
GaAs/AlGaAs heterostructure coated with a film of type-II superconductor jHHj (in 
the latter case the Aharonov-Bohm fiuxes are arranged in a honeycomb lattice, the 
so-called Abrikosov lattice). 

As for the spectral properties of the operator //min, they have been investigated 
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recently in detail by Melgaard, Ouliabaz and Rozenblum In particular, these 
authors proved with the help of results from [Slj and [91J that -ffmin has no zero modes 
at least for periodic lattices of Aharonov-Bohm solenoids, and therefore it differs from 
H^^^ and -f^^ax generic values of magnetic fluxes (and even it is not unitarily 
equivalent to these operators). Let us note that it is possible to extend this result to 
a chain of Aharonov-Bohm solenoids. 



3. The Pauli operator with a singular magnetic field 



In what follows we consider the motion of an electron with the gyromagnetic ratio 
g = 2, therefore 



H 



2m* 



ch 



idx + —A^ ) +[idy + —Ay ) -a^B 



ch 



Let us denote for simplicity 



4A = a, ^B = h. 

ch ch 



(3) 



(4) 



so that dxtty — dya^ = b. In order to employ the dimensionless units we shall consider 
the operator 

H ^ H{^) = -^H{A) . (5) 



h^ 

Introducing a quantum of the magnetic flux, 

27ich 
= 



we also have 



a = —A, b = —B , 
$0 ' $0 ' 



(6) 
(7) 



H = H{a) = {id, + a,f + {idy + ayf - a,b . (8) 

The operator H (and respectively the operator H) decomposes in a sum of two scalar 
operators, 

= H^isi) = [id, + a^f + [idy + ayf ^ 6 , (9) 

(respectively H^{A) = H"^) acting in L^(R^). We admit the vector potential a to have 
singular points, more precisely, we assume that 

a, , a, e LlM^) n C°"(M' \ fi) (10) 

where f2 is a discrete subset (possibly finite or empty) in M? . Consequently, the mag- 
netic field b = dxtty — dytt, is, in general, a distribution in whose singular support 
is contained in Q. Expressions (P) and © represent symmetric operators with the 
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domain C^(M^ \ Q); these operators will be denoted H^{A, Q) and H^{a., Q), respec- 
tively. If the singular support of B coincides with Q (in this case Q is determined by 
the vector potential A) we shall simply write H^{A) and if^(a). 

It is important to note that also in the case when 6 is a distribution the operator 
if^(a) depends, up to unitary equivalence, only on b. More precisely, we have the 
following proposition. 

Theorem 3.1 (gauge invariance of the operator H^{a.)). Let a and a be vector 
potentials with the same magnetic field b (i.e, a, a G Lj'jjj.(M^; M^) fl C°°(R^ \ fi;M^) 
and d^tty — dya^ = d^dy — dyd^ = b in the sense of distributions). Then the op- 
erators H^{a.,Q) and H^{sl,Q) are unitarily equivalent. In more detail, there exists 
a real-valued function f belonging to C°°(]R^ \ fl) such that a = a + grad/, and 
H^{a.,Q) = W~^H^{sl,Q)W where W is the unitary operator acting via multiplica- 
tion by the function exp(— i/). 

Of course, this theorem is well known in the case when the field 6 is a function 
(not a distribution). In the case when 6 is a distribution the theorem is a consequence 
of the following lemma whose elementary proof was communicated to us by K. V. 
Pankrashkin. 

Lemma 3.2. Assume that a G Li^qj,(R^) nC°°(R^ and the equality dxdy — dyUx = 
holds true in in the sense of distributions. Let u E Q and let Q be a rectangle 
containing uo but no other points from Vt. Then 



ax dx + aydy = 0. (11) 

dQ 

Proof. Let us choose functions ip, ip E C^(]R^) so that u ^ supp*^, (p{x,y) = 1 in 
some neighborhood of the boundary dQ, ip{x, y) = (p{x, y) on \ Q and ip{x, y) = I 
on Q. Using the Green formula we obtain 

J axdx + aydy = J ipax dx + ipay dy = J J {dx{fay) — dy{ipax)) dxdy 

dQ dQ Q 

= J J LP {dxtty - dyax) dxdy + j j (oydx^p - axdyip) dxdy 
Qnsupp(i^) Q 

{aydx{ip - ^) - axdy{ip - ^j)) dxdy = . 



Here we have used the fact that the expression dxdy—dyax represents a smooth function 
on \ i7 which necessarily vanishes on this domain. □ 

Proof of Theorem \3.1\ From Lemma we derive in a standard manner that if dxay — 
dyQx = on R^ in the sense of distributions then there exists a real-valued function 
/ G C°°(R^ \ Vt) such that a = grad/ on R^ \ VL. Consequently, if a and a obey 
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the assumptions of the theorem then for some function / G C°^(]R^ \ Q) we have 
a = a + grad/. Let us denote by W the operator acting via multiphcation by the 
function exp{—if). Clearly, is a well defined unitary operator in L^(R^). More- 
over, W leaves invariant the subspace C^(]R^ \ A simple computation shows that 
W~^H^{a, Q)W = i/'^(a, Q). Hence the operators H^{a, Q) and i/'^(a, Q) are unitar- 
ily equivalent. □ 

Remark 3.3. Clearly, if a = grad / in the sense of distributions then dxdy — dya^ = 
in the same sense. 

Remark 3.4. A proposition analogous to that of Theorem IH.ll is also valid for the 
operator H'^^A, fi). Namely, if dxAy — dyA^ = dyA^ — dxAy = B then A = A + grad / 
and H^{A,n) = W-^H^{A,n)W where W = exp{-{ie/ch)f). 

Owing to the gauge invariance it is possible to require the vector potential A to 
have some additional properties. For example, the vector potential A can be frequently 
chosen so that it fulfills the Lorentz gauge condition 

divA = 0. (12) 



4. Basic examples 

In this section we recall several basic examples of magnetic fields fulfilling condition 
()10j). At the same time, we introduce the necessary notation. The majority of results 
presented in the current paper concern Examples 5, 6 and 7. In what follows it will be 
convenient to identify the Euclidean plane with the complex plane C and to work 
with the complex coordinates z = x + iy and z = x — iy. 

Example 1. The homogeneous field 

In this case B = const by definition and one can set 



, B B 



(the symmetric gauge). In the complex coordinates we have 

5 ^ _ , B ^ _ 
Ax = — Im z , Ay = — Kez . 
2 ' ^ 2 

In this example b = 2ttC, where ^ is the number of magnetic flux quanta through a unit 
area in (the flux density). The Lorentz gauge condition (fT^ is obviously fulfilled. 

Example 2. The magnetic field of an Aharonov— Bohm solenoid 

Here B[r) = $(5(r) where $ is the magnetic flux through the solenoid. In this case 
one can set 

A $ y $ a; 
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Equivalently, 

= 9 Im - , tty = 9 Re - , 

z z 

where ^ = $/$° is the number of magnetic flux quanta through the Aharonov-Bohm 
solenoid. Actually, it is well known that 

A \n{\z\) = 2'k6{z). 

In the local coordinates we have 

B = —A - —A = — (—Re--— Im-] = — (— + — ]\n\z\=^6(z) 
dx ^ dy ^ 271- \dx z dy z J 271 ydx"^ dy"^ J 

The vector potential a can be also written as 

a = 6' sgradln |z| . (13) 
Here and everywhere in what follows sgrad stands for the symplectic gradient, 

sgrad/= . (14) 



dy dx 

Hence h = 2ti95{z). The equality diva = trivially follows from |T^. 



Example 3. An arbitrary system of Aharonov— Bohm solenoids 

Let now i7 be a discrete subset of the plane M? and let ($(^)ajgn be an arbitrary family 
of real numbers with indices from VL. We shall consider a system of Aharonov-Bohm 
fluxes intersecting the plane in the points from the set VL and perpendicular to the 
plane. The number equals the flux in the solenoid passing through the point uj 
Then 

h='^^,B = 2T.Y.^J{z-u^) 

where, of course, 9^^ = $^^7$° is the number of magnetic flux quanta through the 
solenoid uj. For a vector potential a fulfilling the Landau gauge condition one can 
choose a meromorphic function M{z) with the following properties: 

1) M{z) has simple poles only, 

2) the set of poles of M{z) coincides with Q, 

3) the residue of M{z) at the point u equals 9^. 

According to the Mittag-Leffler theorem such a function always exists. The computa- 
tions carried out in Example 2 (jointly with the Cauchy-Riemann conditions) show 
that one can set 

ax{z,z) = lmM{z) , ay{z, z) = Re M{z) . 

The operator H^{sl) will be also denoted by the symbol H^{Q,Q) where G = 
{9u))uien- The couple (fi, 6) determines the operator H'^^Q, 6) unambiguously up to 
unitary equivalence. 



10 



Example 4. An arbitrary system of Aharonov— Bohm solenoids with fluxes 
taking a finite number of values 

Separately we consider the case when the number of mutually different fluxes in the 
family ($a;)(^er2 is flnite (equivalently, the family {9^)^^q contains only a flnite num- 
ber of mutually different numbers 6'^). We start from the case when all the involved 
solenoids carry the same flux: 9^^ = 6*, Vu; G Q. In this case we always set 

Here the function W{z) differs from the Weierstrass canonical product Wq{z) related 
to the set Q only by a multiplier exp{g{z)) where g{z) is an entire function. Obviously, 
the set of poles of the function W'{z)/W{z) coincides with fl, all the poles are simple 
and all the residues are equal to 1. Thus one can set 

a = ^sgradln(|iy(2)|). (15) 

Actually, locally we have 
and analogously, 



^Mmz)\) = -lra^'^'^ 



dy " ' W{z) 

In general, let f^i, Qn be mutually disjoint discrete (possibly empty) sets, and 
let 9j, j = 1, . . . , N, he (not necessarily distinct) real numbers. The vector potential 
a is deflned unambiguously, up to gauge equivalence, by the expression 

a = ^ sgradln(|lV,(2)|) = sgradlnf JJ \Wj{z)f' J (16) 
i=i ^ i=i ^ 

where Wj is an entire function having simple zeros only and with the zero set be- 
ing equal to Qj. The function Wj differs from the Weierstrass canonical product re- 
lated to the set Qj only by a multiplier of the form exp{gj{z)) where qj is an arbi- 
trary entire function. An Aharonov-Bohm potential of the form (|T6|) will be called 
a potential of finite type. The operator H^{sl) will be also denoted by the symbols 
H^{n,, . . . , fi^; ^1, . . . , ^;v) or i/±((l^,) ; (%)). 

The most important particular cases of potentials of flnite type are those for which 
the Aharonov-Bohm fleld is invariant with respect to a discrete group A which is 
formed by motions of the Euclidean plane and whose action on Q is co-flnite. 
First of all we shall be interested in the case when the group A is formed by parallel 
translations. Up to isomorphism, there exist just three groups of this type in the plane 
and they are characterized by their rank r (r = 0, 1, 2). 



1. r = 0. In this case A = {0} and the set Q is flnite. 
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2. r = 1. In this case A is isomorphic to Z and has the form A = {kuiQ] /c G Z} where 
ujq is a nonzero vector from . The set Q has the form Q = K + A where K is 
a finite subset of the "elementary strip" F = {x e M^; Q < x-luq < [tuoP} (or, in 
the complex coordinates, F — {z E C; < Rbzloq < |a;oP})- Since each a; e Q 
is uniquely expressible in the form u = k, + X, with k E K and A G A, every A- 
invariant family 6 is unambiguously determined by its subfamily 6^ — {O^jKeK- 

3. r — 2. In this case A is isomorphic to and has the form A = {kiUi + 

k2^2', ki,k2 G Z} where Ui, 002 are linearly independent vectors from M?. The 
set Vl has the form VL = K + K where is a finite subset of the elementary 
cell F = {tiUJi + t2UJ2] < ^1,^2 < !}• We shall assume that the basis a;i, UJ2 is 
positively oriented so that a;i A a;2 = lva.ujiUJ2 > 0. This expression is nothing 
but the area S — S^ol the elementary cell F of the lattice A. 

We shall discuss each of these cases separately. 

Example 5. A finite number of Aharonov Bohm solenoids 
Let A = {0}. In this case the set VL is finite, VL = {a;i, . . . , and 

n 

b^2nJ20,S{z-u;j). 

i=i 

As a rule, the vector potential in this case will be chosen in the form 

n 

a = 9j sgradln(|2; — Uj\) . 
i=i 

The operator i?^(a) will be also denoted by H^{u;i, . . . , Un, Oi, . . . , On). 



Example 6. A chain of Ahcironov— Bohm solenoids 



Assume now that the rank of A equals 1. Firstly we consider the case when K contains 
only one element. Without loss of generality we assume that K — {0}. Then A, 
= 6 for all a;, and 



wu(z)-^ n (i-£:)''"""=^n(i 

fcGZ.fe^^O ^ k=l ^ 



Therefore one can set 
Consequently, 



kW^ 



1 

— sm 

TT 



W{z) = sin( — 



a = 9 sgrad In 



sm 



nz 
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which means that 



Ox = — Im ctg — , a„ = — Re ctg — . 
Generally, Q, — K -\- K with an arbitrary finite subset K C F, and we have 



Then the vector potential reads 

' TT 

In I sml 

KeK 



a = sgrad 9^ In ( sin I 



U — K 



Example 7. A lattice of Ahcironov— Bohm solenoids 

Assume now that the rank of A equals 2 which means that A is a two-dimensional 
lattice. Again, we shall start from the case when K = {0}, hence = A. In this case 
Wa{z) coincides with the Weierstrass cr-function of the lattice A, 



n (l--)axp(i + 4 

At the same time. 



a[z;uji,uj2) = a[z) = z 

u>Gn\{o} 



a{z) 

is the Weierstrass ^-function of the lattice A. Thus 

a = ^ sgrad HH^)\) = ^ C^^Ci^), ^eC{z)) . 

In the general case Q — K + A with an arbitrary finite subset K C F. Then the 
magnetic field takes the form 

KeK AeA 

One can set 

a = sgrad hi{\a{z — k,)\) . 

KeK 

Remark 4.1. In all the examples with an Aharonov-Bohm potential of finite type, 
and also in the case of a homogeneous magnetic field, there exists a function (fi{x, y) 
such that a = sgrad (f or, equivalently, 

Aip^b (17) 
Namely, one can respectively set in Examples 1, 2, 3, 5, 6, 7: 
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ip{x,y) = -e \n{x' + y'') = eH\z\) 



ip{x,y) = In \W,{z)\ = In J] I^^ WI' 



0=1 



¥^(2;,?/) = X^^i ln(|2 -c^jl) 



/ 7r(z — k) 
sin ' 



n 



sin 



= ^ ln(|o-(2; - k;uji,uj2)\) = ln( JJ |o-(2; - 



Let us note that in the general case when is a A-periodic continuous field the 
solution of the equation (|T7j) is expressible in the form 



(18) 



^{z) = — II H\a{z-z')\)b{z')dx'dy'. 



where F is an elementary cell of the lattice A Actually, we have already seen that 

Aln(|cT(z)|) = 27r ^ - A) . 



Therefore a formal computation yields 

Aip{z) =J2 [[ ^^^'^ -z' -\) dx'dy' . 



AeA 



(19) 



For every z G C there exists a unique Aq G A such that 2; G F + Aq, i.e., 2; — Aq G F. 
Then the summands in fll9|l with A 7^ Aq vanish and we have 

A^{z) = jj h{z') 5{z' -{z- Ao)) dxdy = b{z - Aq) = b{z) . 

F 

In the case of a lattice formed by Aharonov-Bohm solenoids formula |TH|) makes still 
sense and it again yields 

viz) = e ln{\aiz)\) . 

Let us note that the Lorentz gauge condition ()12|) follows from the equality a = 
sgrad (p. 
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In the sequel, the main results will be derived for Hamiltonians corresponding to 
three types of systems of Aharonov-Bohm solenoids. Namely, the set fl formed by the 
intersection points of solenoids with the plane may be 1) a finite set, 2) a chain or a 
finite union of chains, 3) a lattice or a finite union of lattices. These systems will be 
called regular. 

5. A rigorous definition of the Pauli operator as a 
self-adjoint operator 

Let us return to the symmetric operators = H^{sl, Q) defined in Q while assuming 
that condition (jTHI) is satisfied. Let us introduce the momentum operators 

= P.^{sL,n) = -idx - , Py = Py{a,n) = -idy - ay . (20) 

In virtue of (fTUI) these operators can be considered as symmetric operators in L^(M^) 
with the domain C^(M^ \ fl). Following Aharonov-Casher we define the operators 

T^ = T±{8i,n) = P,±iPy, (21) 

or T+ = —2idz — A{z, z),T^ = —2idz — A{z, z) where A = ax + iay. Then the following 
equalities hold true on C^iR^ \ Vt): 

T+T_ = H , T_T+ = H+ . (22) 
By a straightforward computation one can verify a simple but important lemma. 
Lemma 5.1. The commutation relations 

[Px,Py]=tb, [T_,T+] = -26, (23) 

are valid on C^(M^ \ ^)- -^^ particular, z/ supp B G Q (including the case when B 
corresponds to a system of Aharonov-Bohm solenoids ) then the operators Px and Py 
(respectively T+ and T_) commute on the domain C^i^ \ Vt). 

iFiom. the obvious inclusions 

Tl D (24) 

we immediately deduce that the operators T± are closable and therefore the self-adjoint 
operators 

^min = ^min(a,^^)=TiT± (25) 

are well defined (see, e.g., [HZl Theorem X.25]). The associated quadratic forms h^^^ 
are closures of positive forms defined on C^(R^ \ Q) by the expressions 

(T±(^|T±^), (26) 

respectively. 
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On the other side, let us consider a quadratic form defined on C^(]R^ \ Q) by the 
relation 

S^iip, ^) = { PMPx^ ) + ( Py^lPyi^ ) T ( bifl^lJ ) . (27) 

By a straightforward computation using relation one can show the following 
lemma. 

Lemma 5.2. The quadratic forms /i^j^ and coincide on C^(]R^ \ il). 

In particular, if the support of B is contained in Q then the quadratic forms /ij^j^ 
and coincide on C^(M^ \ Q) and therefore they are necessarily equal. 

Corollary 5.3. If B is a distribution with a support contained in Q then the operators 
-^min ^'^^ -^min coincidc. In particular, for the vector potential a of a system of the 
Aharonov-Bohm vortices supported on a set Q cM."^ these operators coincide with the 
Friedrichs extension of the symmetric operator defined on C^(]R^\fi) by the differential 
expression 

i/o(a) = (id., + a.,f + {idy + ayf . (28) 

In view of Lemma (5.21 we shall sometimes simply write -f^min instead of -f^min- The 
operator -fTmin has been investigated in detail in [89. . 

Jointly with the operator iJmin(a, f2) let us consider the operators 

^max = ^max(a,^^)=T^T; (29) 

with the associated quadratic forms defined on V{H^^^) by the expressions 

hLM^) = {T;^\T;^), (30) 

respectively. 

The definitions of H^^^{sl, fi) and H^-^^{ai, fi) in principle depend on the choice of 
the discrete set Q. If Q coincides with the singular support of b, however, we shall 
simply write, similarly as in Section El H^-^^{ai) and H^^^{sl) since in that case the 
vector potential a determines Q unambiguously. 

If the field B is sufficiently regular and = then the operator H^-^^ coincides 
with the operator -ffj^ax This is not true, however, for operators with Aharonov- 
Bohm fiuxes (see jHj, [Hnj)- Since in this case -ff^in defined by expression 
and is independent of spin, this operator is the Schrodinger operator of a spinless 
particle in the presence of the Aharonov-Bohm fiuxes (or the Schrodinger operator of 
a particle with spin when interaction of the spin with the field can be neglected). On 
the other hand, -f^min defined by expression ©, they depend on the spin and may 
be considered as the Pauli operators for an electron with the gyromagnetic ratio g = 2. 
Below we are interesting in the properties of ground states of the operator -f^^ax- 

For the analysis of operators -f^^ax following description of the operators 
will be useful. Namely, owing to condition ^TU^ the differential operators —id^ — and 
—idy — ay are well defined on the space of distributions V'{M? \ Vt). Consequently, the 
operators T± defined on C~(M2 \ ^) 

can be naturally extended to linear mappings T± 
defined on D'(]R2\fi). Using the fact that L'^{M?) is naturally embedded into V'{M?\n) 
we get the following lemma. 
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Lemma 5.4. The operator is a restriction ofT^ to the domain 

{f e L\R')- f^f e L\R')} . 

Using this observation we can prove the following lemma. 

Lemma 5.5. Let C he the operator of complex conjugation, Cf = f. Then CH^.^^{sl, Q) 
/7^,,(-a, n)C and CH^Je., fi) = H^J-^, n)C . 

Corollary 5.6. The operators H^^^{sl,Q) and H^^^{—a,Q) have the same spectra. 
In particular, they have the same eigenvalues with equal multiplicities. An analogous 
proposition holds true for the couple of operators H^^^{sl, Q) and a, Q). 



6. Elimination of Aharonov— Bohm solenoids with 
integer fluxes 

In this section we consider a vector potential a of the form 

a = a + a.AB 

where a^^ is a vector potential corresponding to a system of Aharonov-Bohm solenoids 
intersecting the plane in the points of fl. We describe here briefly the "gauge-periodicity" 
of the operators with the vector potential a; details can be found e.g. in jT^l IHOl EH E21 • 
First we shall assume that the considered solenoids carry equal fluxes of the value 
9ab- In this case we set aab = Qab sgrad ln(|W^(z)|) (cf. Example 4 in SectionE]). Let 
9ab be an integer. Then the function 

g{z,z) = exp(^9AB ln|^^|^^^ = exp{i9AB aTg{W{z))) , 

is well defined and continuous in the domain C \ Clearly, \g{z, z)\ = 1, E C\Q, 
and, moreover, g G C°°(]R^ \ Q). 

Lemma 6.1. If 9ab is an integer then the following relations hold true 

g-'P,{k, n)g = P,(a, ^]) , g~'Py{k, n)g = P,(a, ^]) . 
Proof. It suffices to show that 

-i grad g = g SLAB ■ (31) 

Actually, we have 

o-x'iwi\) ^ {l,^^E){^^m^))-l{Hwiz)) + Hwi^)) 
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Analogously, 



d f W{z) \ W'{z) _i 



Relation obviously follows from these equalities. □ 

Assume now that slab is a vector potential corresponding to an Aharonov-Bohm 
field of finite type whose singular support coincides with i7 = fii U . . . U and with 
an array of fiuxes denoted by 9 = {Oj)i<j<n- Let {jnj)i<j<n be an arbitrary array of 
integers and let a^^ be another Aharonov-Bohm potential of finite type defined by 
the same array of sets {Vtj) but with the fiuxes 9j = 6j + rrij. 

Theorem 6.2. Assume that a, slab one? a^^ have the same meaning as described 
above. Then H'^^^{ai + slab, ^) (respectively, H^^^{8i + slab,^)) is unitarily equivalent 
to the operator iy^jj^(a + a^s, Vl) (respectively, H^^^{sl + slab, ^) )■ 

Proof. Let Ti±, T2± be the operators corresponding to the vector potentials a + slab 
and a + slab, respectively, as described in Sectional By construction, 

V{T,±)=V{T2±) = C^{R'\n). 

Applying repeatedly Lemma f6. II one can show that there exists a unitary operator U 
such that 

f/(Co~(M2 \n))= c^im" \ n) 

and 

U-'T2±U = Ti±. 
From the unitarity of U it follows that 

U-^T2±U = Ti± and U-^T;JJ = T*±. 

Consequently, 

U-'H^,^{SL + SLAB, m = U-'T;^T2±U = T*^T,± = if± ,(a + a^^, Q) 

and 

U-'Ht.{sL + SLAB, m = U-'T2^T;^U = T,^T*^ = Ht^isL + slab, n). 

This shows the theorem. □ 

Corollary 6.3. If all fluxes 6j are integers then the operator H^^^{8i+ slab , ^) (respec- 
tively, H^^^{sl+slab, ^) ) is unitarily equivalent to the operator H^-^^{sl, fi) (respectively. 

Let us formulate separately two most important cases of this corollary. The first 
one is based on the fact that the both operators iJj^i^(0, fi) and H^^^{0, ^l) do not 
depend on the choice of the discrete set Q and coincide with the Laplace operator 
-A. 
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Corollary 6.4. Let Q, be a discrete set which is invariant with respect to a co-finite 
action of a lattice A of rank r, < r < 2. Assume that a = and slab is a vector 
potential corresponding to a system of Aharonov-Bohm solenoids supported on the set 
Vl and such that all fluxes are integers. Then each of the operators ifj^jjj(ayiB, Vl) and 
^m&A^ABi^) is unitarily equivalent to the Laplace operator —A. 

Proof. Since the action is co-finite we are again in the situation when Q sphts into a 
finite union = i^iU. . .Uf2„. Hence one can apply Corollary 16. 31 The unitary operator 
induced by multiplication with the function g acts locally in the form sense jHSI and 
therefore each of the operators H^^^^(0, fi) and H^^^{0, f2) is a point perturbation of 
—A supported on the set fl. The perturbed operator is clearly positive and local in 
the form sense [HI]. On the other hand, every nontrivial point perturbation in the 
two-dimensional case is known to have a strictly negative infimum of the quadratic 
form over unit vectors [3^] • D 

Since the minimum of spectrum in the case of a periodic point perturbation of the 
Landau operator is strictly smaller than the minimum of spectrum of the unperturbed 
operator |nH| the following corollary is also true. 

Corollary 6.5. Let a be a vector potential of a nonzero homogeneous magnetic field 
and assume again that the discrete set Q is invariant with respect to a co-finite action 
of a lattice A of rank r, < r < 2. Then for b > 0, each of the operators H^^^{a + 
SiAB,^) and H^^^{a + a^B,^) is unitarily equivalent to the Landau operator _ff+(a). 
For b < 0, an analogous statement is true for the operators i/~i^(a + slab,^) and 

To simplify the discussion to follow we shall assume once for all that an appropriate 
gauge transformation has been applied so that the values of all involved Aharonov- 
Bohm fluxes belong to the interval [0, 1[. If there are some zero values then VL is 
strictly larger then the singular support of b. As shown by Corollaries 16.41 and 16.51 
the zero values can be eliminated in some particular cases. We shall proceed in our 
simplifications even further. If not said otherwise, we assume everywhere in what 
follows that the values of Aharonov-Bohm fluxes belong to the interval ]0,1[ and, 
consequently, the singular support of b coincides with VL. 



7. The ground states (zero modes) of the Pauh op- 
erator 

It follows immediately from the definition of the operators -ffmax -ffj^m ^^^^ they 
are nonnegative. Consequently, if the equation 

H^J^ = (32) 



or the equation 

^max^ = (33) 
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has a solution in L^(R^) then this solution ip± (called zero mode) is a ground state 
of the corresponding operator. Since the equality H^-^^ijj = implies { H^^^ijj\ijj) = 0, 
i.e., the equality ||T±'?/'||^ = 0, equation (jH^ is equivalent to the equality 

T±^ = 0. (34) 

Analogously, equation (jH!^ is equivalent to the equahty 

T^7/> = 0, (35) 

or, this is the same, to the condition 

f±^ = 0, ^gL2(M^). (36) 

Suppose that the vector potential a was chosen to have the form a = sgrad if 
where (f satisfies the equation Ay? = 6 in the sense of distributions. We shall seek a 
solution of equation in the form 

^±(a;, y) = exp(=Fv2(a;, y))f{x, y) = exp{T^{z, z))f{z, z) , (37) 

where / has to be chosen so that '\l}± G L^(R^) [the Aharonov-Casher ansatz). In the 
space of distributions ©'(M^ \ Vt) we have 

= exp(-v5) [[i{d^^p - ay)f + {-dyip - a^)f) - i(9x./ + idyf)) 

= -2«exp(-v9) — 

and 

= exp(v?) {{i(-dxip + ay)f + {-dyip - a^)f) - i{d^f - idyf)) 

= -2iexp{(f) — . 

/^From here we deduce that the relation 

ifl>+ = 0, ^+gL2(M2), (38) 
is equivalent to the condition 

^ = 0{zeC\n), exp(-(^)/ G L2(M2) . (39) 
Analogously, the relation 

/7^,,^_ = 0, ^_gL2(m2), (40) 
is equivalent to the condition 

{zeC\n), exp{ip)f e L'^iR'^). (41) 



dz 

This shows the following theorem due to Aharonov and Casher [Hj 



20 



Theorem 7.1. Assume that a vector potential a is expressed in the form a = sgrad (p 
where ip satisfies the equation A(p = b in the sense of distributions. Then solutions 
of the equation H^^^ip = in L^(]R^) are exactly those functions from L^(M^) which 
have the form ip-\-{z, z) = exp{—ip{z, z))f{z) where f is a holomorphic function in the 
domain C\Q. 

Similarly, solutions of the equation H'^^ip = in L^(M^) are exactly those func- 
tions from L^(R^) which have the form ip^^z.z) = exp{{p{z,z))f{z) where f is a 
holomorphic function in the domain C \ ^2. 

Let us point out an interesting consequence of the theorem. 

Proposition 7.2. Assume that the both operators H^^^ and H^^^ have zero modes. 
Then they are distinct. In particular, the set C^(C \ Q) is not a core for at least one 
of them. 

Proof. Let be a zero mode of H^^^. Suppose that this operator coincides with 
-^max- Then ip is a. zero mode for H~.^^ as well. Using notation of Theorem 17.11 we 
have ijj = exp{—ip)f = exp{ip)g where / is holomorphic in the domain C \ and g 
is antiholomorphic in the same domain. Since ip is real it holds true that {ipl"^ = fg. 
Taking into account that g is holomorphic the last equality implies that fg is a constant 
function and hence the same is true for I?/' p. Since G L^(]R^) it follows that ip = 0, 
a contradiction. □ 



8. Zero modes of the operators H^^^ with Aharonov— 
Bohm potential of finite type 

8.1. Formulation of the problem 

In this section we shall study ground states of the operator H^^^{a, Q) for an Aharonov- 
Bohm potential a of finite type determined by mutually disjoint discrete sets f2i , . . . , f2„ 
such that = fii U . . . U Qn, and by fluxes (not necessarily distinct) 9i , . . . ,9n (cf. 
Example 4 from Section Recall that we assume that < < 1, for all j. 

We can rephrase the formulation of the problem. Namely, according to Theo- 
rems 16.21 and 17.11 we have to study square integrability of a function ijj having the 
form 

n 

^{z,z) = f{z)l[\W,{z)\-'^+"^^ (42) 

i=i 

where the numbers rrij are integers, f{z) is holomorphic or antiholomorphic in the 
domain C \ ^2, and the functions Wj determine the potential a according to formula 

(CHI)- 

In this section the following lemma will be useful. 
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Lemma 8.1. Assume that a function f{z) is expressible in an annulus ri < \z\ < r2 
as a Laurent series 

CO 

n=— oo 

Then for ri < r < r2 and arbitrary n & Z it holds true that 

1/(^)1 |dz|>27r|a„|r'^+^ (43) 



\z\=r 

In particular, if ri = and n > —2 then 



11 \f{z)\dxdy> 



27r|an| 
n + 2 



,n+2 



\z\<r 



if ri = and a„ 7^ for some n < —2 then 

\f{z) \ dxdy = 00 



\z\<r 



if r2 = 00 and a„ 7^ for some n > —2 then 



\f{z) \ dxdy = 00 . 



\z\>r 



Proof. The proof immediately follows from the simple estimate 



27r 



2n 



1/(^)1 |d^|=r J \f{re'n\dv>r 

\z\=r 

An application of Lemma f8 . II yields the following auxiliary result. 
Lemma 8.2. Assume that in it holds rrij = 0, for all j , i.e., 



ij{zrz) = f{z)\{\w, 



(44) 



(45) 



(46) 



27r|a„|r"+^ □ 



where —l<(3j<Q, for each j , and f{z) is holomorphic in the domain C\Q. If 
ip G L^(]R^) then f is an entire function. 

Proof. Assume that u & and r > is sufficiently small so that D{u, r) (IQ = {ui} 
where D{uj,r) is the disc with radius r centered at the point u. Since Wj^uj) 7^ 0, for 
j 7^ k, and u; is a simple zero of Wk{z) one can assume that r is small enough so that 
it holds 

n 

n > c> 0, 
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for < |z — < r. 

Now one can show that u cannot be a pole nor an essential singularity of /. 
Otherwise uj would be a pole of even order or an essential singularity of In any 
case Lemma f8. II implies that 

p{z,z)\'^dxdy>c'^JJ \f{z)\'^ dxdy = oo , 

D{uj,r) D{uj,r) 

a contradiction. □ 



8.2. Finite number of Aharonov— Bohm fluxes 

We start from the simplest case, i.e., from the Hamiltonian H^^^{Q, G) corresponding 
to a finite number of Aharonov-Bohm fluxes (Example 5 from Section 0]). Then Q = 
{ui, . . . , Un} is a finite set, 9 = (^i, . . . , 9n), < 9j < 1. In this case zero modes may 
occur under suitable assumptions on fluxes 6j. More precisely, the following theorem 
is true [Hj . 

Theorem 8.3. A sufficient and necessary condition for the operators H^^^{Q, 0) and 
-^max(^' 6) to have zero modes is 

n 

in the former case, and 

n 

E < n - 1 (48) 

in the latter case. 

Proof. Let us start from the operator -ff^ax- have to find a nonzero function / 
which is holomorphic in the domain C\Q and such that the function 

n 

ij{z,z) = f{z)l[\z~u,\-'^ (49) 
i=i 

is square integrable. Suppose that condition (jTfj) is satisfied. Taking for / a constant 
function it is easy to verify that in that case we get a square integrable function ip. 

Conversely, assume that is square integrable but condition (jTfj) is false. Then 
from ()49p one easily deduces that / cannot be a nonzero constant. Furthermore, from 
the equality 

f{z) =iIj{z,z) Yl \z-u;jf^ , 
we find that there exists a constant ci > such that 



1/(^)1 <ci(l + |^|)|^(;^,z)| onC. 
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Consequently, if r > 1 then 

J J \fiz)\'dxdy<C2r' 

\z\<r 

where 

C2 = 4c^ // \^jj{z^ z)]"^ dxdy . 



From inequahty ()43|1 in Lemma l8.ll it follows that f{z) is a constant function. This 
contradiction proves the theorem in the case of the operator H^^^. To prove the 
theorem in the case of the operator H~^^ one can either modify the above argument 
or to apply CoroUarv 15.61 □ 

8.3. A chain of Aharonov— Bohm fluxes 

Here we show that the Hamiltonians -ffj^ax corresponding to a finite union of chains 
of Aharonov-Bohm fluxes have infinitely many zero modes. In this subsection we use 
the notation from Example 6 in Section 0] 

The proof uses the following elementary estimate. Since for z = x + iy, x, y & M., 
we have | sin(2;)p = ch^{y) — cos^(x) = sh^(?/) + sin^(a:), it holds true that | sh(?/)| < 
I sin(^)| < ch.{y). Hence 

< |sin(z)| <el^l. (50) 

Theorem 8.4. Let a uniformly discrete set Q be expressible as a disjoint union of a 
finite number of chains Qi , . . . ,Qn, o-nd let the chain VLj = Kj + Aj carry Aharonov- 
Bohm fluxes (6'k)kga'j (j = I , . . . , n, < 9^, < 1). Then the Hamiltonians H^^^{Q) 
have infinitely degenerate zero modes. 

Proof. In the proof we shall consider only the operator i^^ax- Using Lemma lA.ll we 
may assume that each chain Qj is contained in a line Lj and that it holds Lj ^ Lk for 
j 7^ k. Then the Bravais lattice of the chain Qj has the form Aj = {kuj] k G Z}, 
Uj G C, ujj 7^ 0. Without loss of generality we can suppose that ui > and ki = 0. 
Hence Li = M and Lj = ujj'R + kj (j = 2, . . . , n) where Kj is a fixed element from Kj. 

For each line Lj we shall construct a strip Pj with border lines parallel to Lj and 
containing Lj in its interior. Furthermore, let Q be a sufficiently large disk centered 
at such that outside Q the strips Pj do not intersect each other. 

It suffices to show that there exists an infinite number of linearly independent 
entire functions f{z) for which the function 

z) = f{z)9i{z, z)- gn{z, z), 

with 

(51) 



9j{z,z)= TT sin(^^^^ — —] 

.at- V J 
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is square integrable. We shall show that this condition is satisfied for any function 

sm(az) 



where 



0<a<^:= — V^^. (52) 

The verification follows from a series of claims. 

(A) ^ G L2(Q), 

(B) each function gj is bounded outside the strip Pj, 

(C) the function gi{z, z) sin(a;2;) is hounded outside the strip Pi. 

Claim (A) follows from the fact that / is bounded on Q and that the functions 
Qj have square integrable singularities. Moreover, only finitely many singularities are 
contained in Q. Claims (B) and (C) are consequences of the inequalities in (jSUj) and 
condition (|^. 

To complete the proof it remains to show that 

(D) V^GL2(P,\g),Vj = l,...,n, 

(E) V'eL2(M2\(PiU...UP„)). 

To show Claim (D) notice that (B) and (C) imply the estimate 

\^{z,z)\ < ^ \gj{z,z)\, 

valid on Pj \ Q with some constant Cj > 0, and that the function gj{z,z) is periodic 
along the line Lj. For j = 1 one uses also that sin(az) is bounded on the strip Pi. 
To show Claim (E) let us point out that the inequahty 



(z, z)\ < c 



simaz] 



\9i[z,z) 



holds true on C \ (P2 U . . . U P„) with some constant c' > 0, as it follows from (B). 
From inequalities (|^ one derives that 

on C \ (Pi U P2 U . . . U P„). Finally, condition §^ implies that ^jj e L^iR"^ \{PilJ . . .U 
Pn)). □ 

Under more restrictive conditions on the fiuxes the assumption on the uniform 
discreteness can be dropped. Since every chain is a union of one-atom chains we can 
confine ourselves to such chains. Moreover, it is clear that a union of chains need not 
be a uniformly discrete set only in the case when among the chains in question there 
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are at least two contained in the same line. Consequently, it suffices to analyze the 
case when the chains are contained in a single line, say, in the real axis R. Suppose 
that Qj = Hj + Aj where = {ujk; k G Z}, Hj G M, ujj > (j = 1, . . . ,?t,), are 
mutually disjoint one-atom chains and 6'i, . . . , 6'„ (0 < 9j < 1) are the corresponding 
Aharonov-Bohm fluxes. 



Theorem 8.5. Assume that all chains Qj, j 



n, are contained in 



Then 



the Hamiltonian -f^^ax = -^max(^i5 • • • 5 ^n] di, ■ ■ ■ , dn) has an infinitely degenerate zero 
mode if one of the following conditions is satisfied: 

(i)ei + ... + en<i, 

(llj h . . . H > h . . . H : , 



(Hi) n = 2. 

Proof, (i) In this case one can choose numbers pj > 1 so that pjBj < 1, Vj = 1, 
and J2pJ^ = 1 (s-g-, Pj^ = dj/ (^^i + . . . + 9n))- Let us consider the functions 



9jiz,z) 



sm{aj{z-Kj)) 



Z Kj 



sm 



7r(z — Kj) 



Uj 



where 



vr 



< < min — Qj . 



3 iOi 



Set ip = Qi - ■ ■ Qn- It suffices to show that ip is square integrable. From the Jensen's 
inequality it follows that 



< 



1^1 



2pi 



Pi 



+ . . . + 



\97 



\2p„ 



Pn 



(53) 



Recalling that Pj9j < I, pj > 1, and repeating the considerations from the proof 
of Theorem 18.41 one can show that each summand on the RHS of formula ()53|) is 
integrable. 

Let us now discuss condition (ii). One can assume that minjUj = ui and ki = 0. 
We shall consider a function ip of the form 

smjaz) 

^[z, z) = gi{z, z)--- gn{z, z) 



where 



9i{z,z) 



sm — 



(54) 
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for j = 2, . . . , n, and a obeys the condition 

0<a<e:= 7r> ^. 56 

3 ■' 

Let T be a strip parallel to the real line M and containing M in its interior. One again 
concludes from (|Kn|) that outside the strip T it holds true that 

Furthermore, inside T we can use the estimate 

\i){z,z)\ < C2 



sm a2; 



Therefore one can use a similar reasoning as in the proof of Theorem 18.41 to show that 

Finally, let us discuss condition (iii). If ui = 002 then one can refer to Theorem 18.41 
In the opposite case we shall assume that Ui < 002- Oi + O2 < 1 then we apply 
condition (i) from the theorem. If not then we have 

gl ^ ^2 ^ ^1 + ^2 ^ 1 _ 1 ^ 1 1 



uji UO2 0J2 0J2 u!i 0J2 mm UJj 

and we can apply condition (ii). □ 

According to Lemma 15.51 we can reformulate the result for the operator -ff~ax 
follows. 

Theorem 8.6. Assume that all chains Vtj, j = l,...,n, are contained in R. Then 
the Hamiltonian -f^^ax — -^max(^i; • • • > ^n! ^1, • • • > 6'„) has an infinitely degenerate zero 
mode if one of the following conditions is satisfied: 

(i) 6^ + ... + en>n-l, 
(11) h . . . H < 



uJi Un mm Uj 

j 

(iii) n = 2. 



8.4. A lattice of Aharonov— Bohm fluxes 

Let us now consider the Hamiltonian -ffj^ax ^r a lattice of Aharonov-Bohm solenoids 
Q = K + A where A is the Bravais lattice of the crystallographic lattice Q with a 
basis {ti;i,ti;2} (cf. Example 7 in Section HJ. To analyze this case we shall use the 
Weierstrass cx-function a{z) = a{z;u!i,u!2)- Let us introduce, following jHSl EH], the 
modified Weierstrass cr- function a{z), 



(57) 



27 



where 



S = lm{uJiLJ2] 



and ({z) is the Weierstrass ^-function (cf. Appendix [n)) . We shall need the following 
lemma. The number /i, occurring in the formulation of lemma and depending on the 
lattice A, is defined by 

(58) 



Lemma 8.7. Let aj, j = 1,. . . ,n, be real numbers such that < aj < 1, let (3 be 
an arbitrary real number, and let aj, j = 1, . . . ,n, be an arbitrary array of complex 
numbers such that among them there is no couple congruent modulo A. If the condition 



a A 



is satisfied then 



exp(p|^| 



n 



Proof. We shall consider a shifted elementary cell of the lattice A, 

Le = {{ti + e)uji + {t2 + e)uj2; < ti, t2 < 1} , 

where e > is chosen so that the interior of contains exactly one zero for each of 
the functions a{z — aj), and hence exactly one pole of the function l/cr{z — aj). But 
in that case, 

-2a i 



n 



aiz ~ a,- 



dxdy < oo . 



Let p{z, z) be a function defined by the formula 

|cr(2;)|^ = exp(i/2;^ + uz"^ + 2iJLZz)p{z,z) 
From formula (|Un|) we deduce that 



JJ^ \p{z — aj, z — aj) I dxdy < oo 
i=i 



(59) 



(60) 



(61) 



Since the function p{z, z) is A-periodic and |cr(2;) ^ = exp(2/^|2;p)p(z, z) (see Lemma lTTII 
it holds true that 



exp(2/?|2| 



z — ai 



dxdy 



= E //exp(2(/3-/.5^ 
< /^sup{exp(2(/5-/i5^ 



aj)\z\ 



Ylpiz -aj,z- aj) 



dxdy 



aj)\z\ 



z e Ls + X> < oo . □ 



aga 
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Theorem 8.8. Let = K + A be a lattice of Aharonov-Bohm solenoids with an array 
of fluxes = {6k)k&k, < < 1. Then each of the operators H^^-^iVL^Q) has an 
infinitely degenerate zero mode. 

Proof. We shall confine ourselves to the case of the operator -ffmax- Let us consider 
the function 

^{zrz) = f{z)\[\~a{z-K)\-''^. (62) 

According to Lemma f8. 71 ip G L?{E?) if / is an arbitrary polynomial. □ 

Remark 8.9. Owing to Theorem 18.81 we can describe an interesting example related 
to the question of absolutely continuous spectrum for the Pauli operator H^^-^{a) with 
a magnetic field h = dxdy — dya^ which is supposed to be periodic with respect to a 
lattice A = {kiUi + ^2(^2; ^1, ^2 ^ with S = lm{ujiU!2) > 0. If the vector potential 
a is "sufficiently regular" and the flux of the field b through the elementary cell equals 
zero then the spectrum of the operators H^^^{a.) is purely absolutely continuous (see 
j2Sl 1211 I2ni I2H] and others). The same result is true for Schrodinger opera- 
tors with "sufficiently regular" periodic vector potentials in the space L'^{W^) for any 
d > 2. In the case d > 3, N. D. Filonov described in jSH] an example showing that 
the assumptions on the vector potential stated in [2H] and other papers cannot be 
essentially weakened. Theorem 18 . 81 shows that two-dimensional Pauli operators with a 
singular two-periodic magnetic field may have (infinitely degenerate) eigenvalues. In 
more detail, let us take, for example, a set K containing two elements, K = {ki, K2}, 
and suppose that 9 = 9^^ = —9^2 ^]0,1[. Then by Theorem 18.81 the both operators 
-^max(^) have an eigenvalue, namely the number zero. According to Example 7 from 
Section m the corresponding vector potential a reads = 9 lm{({z — ni) —({Z — K2)), 
Qy = 9 Re{({z — Ki)—({z — K2)). Owing to quasi-periodicity of the Weierstrass function 
C{z), ({z + ujj) = ({z) + Tjj where r]j = 2 Q{ujj/2), the vector potential a is A-periodic. 

Now we state an analog of Theorem 18.41 for lattices of Aharonov-Bohm fiuxes. In 
view of the example described in Remark IA.2[ we cannot here repeat the arguments 
from the proof of Theorem 18.41 but the properties of the modified Weierstrass a- 
function simplify matters considerably. 

Theorem 8.10. Let a uniformly discrete set VL be expressible as a disjoint union of a 
finite number of lattices f2i , . . . , VL^, and let the lattice VLj = Kj + Aj carry Aharonov- 
Bohm fluxes {9f^)^(zKj (j = I , . . . , n, < 9^ < 1). Then the Hamiltonians H^^^^^Q) 
have infinitely degenerate zero modes. 

Proof. In the proof we shall consider only the operator -ff^ax- Without loss of generality 
we suppose that each Kj is a singleton: Kj = {nj}, and we shall write 9j instead of 
6*^ . By the hypothesis, there is a sufficiently small disk D centered at such that for 
uji ,uj2 ^ 7^ ^2, the sets D+uJi and D+UJ2 are disjoint. Denote Lj := D + Hj+Aj , 

then for every j there exists Cj > such that \cr{z — < Cj for z ^ Lj. It is clear 

that 

n 11 

n - a,)n < E ( n ^'^) - ■ 

j=l j=l k^j 
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Now we can refer to Lemma f8. 71 □ 

Analogs of Theorem 18.51 and Theorem IS.fil are also valid and can be proved by 
the same method. In more detail, let Qj, j = l,...,n, be mutually disjoint simple 
crystallographic lattices, f2j = kj+Aj where kj G C, = {kiooi^ +k2UJ2'^ ; ki, k2 ^ Z}. 
Furthermore, Sj = lm{ui[^^uj2 '') designates the area of the elementary cell of the Bravais 
lattice Aj. 

Theorem 8.11. The Hamiltonian H^^^ = H^^^{Qi, . . . ,Qn](^i, ■ ■ ■ ,(^n) has an in- 
finitely degenerate zero mode if one of the following conditions is satisfied: 

(t)e^ + ... + en<i, 

(ll)— + ... + — >— + ...+ 



oi bn oi bn mm bj 



1 



(Hi) n = 2 and Si ^ S2 ■ 

Theorem 8.12. The Hamiltonian H^^^ = H~^^{Qi, . . . ,Qn;6i, . . . ,6n) has an in- 
finitely degenerate zero mode if one of the following conditions is satisfied: 

(i) 01 + ... + 9n> n-1, 

y. ^ 0^ 6r) 1 

TT + • • • + TT < 



bi bn mm bj 



(Hi) n = 2 and Si ^ S2 ■ 



8.5. Superposition of a homogeneous magnetic field with a 
field corresponding to Aharonov— Bohm solenoids 

Here we consider a perturbation of a homogeneous magnetic field by the field corre- 
sponding to a system of Aharonov-Bohm solenoids, i.e., we consider a vector potential 
a of the form a = ao + a^^ where ao is the vector potential of a homogeneous magnetic 
field 60 = 27r^o with a flux density ^0 = '^i^i.—y, x), and S-ab is the vector potential 
of a system of Aharonov-Bohm fluxes. We shall suppose that the potential a^^ is of 
finite type. In that case we have a finite family of mutually disjoint discrete subsets 
in the complex plane, fli, . . . , fin, and in each point of the set Qj (j = 1, . . . , n) there 
is a flux of magnitude 9j (0 < 9j < 1) intersecting the plane. 

Suppose for deflniteness that 60 > 0. Then the operator H^^^{ao) has an inflnitely 
degenerate zero mode (the lowest Landau level shifted by the value —bo) while the 
ground state of the operator H^^{a.o) is strictly positive (this is the lowest Landau 
level shifted by the value bo) - Thus the latter operator has no zero mode. Intuitively, the 
results proved below in Theorems 18 . 1 HI and 18 . Ifjl mean that if the set f2 = f2i U . . . U0„ 
has a flnite density then a superposition with the potential a^B does not remove the 
zero mode from the spectrum of the operator H^^^, and a zero mode cannot occur in 
the spectrum of the operator H~^^ provided bo is sufficiently large. Moreover, if this 
set has zero density then the same statement about zero modes of the operators -f^^ax 
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and i?~ax t'^^^ fo'^ ^^?/ &o > 0. In the case when f2 is a lattice, a superposition with 
the potential a.AB does not remove the zero mode from the spectrum of -ff^ax ^^2/ 
6o > but a zero mode may occur in the spectrum of H~^^ for particular values of 
fluxes 6j. An attentive reader can effortlessly guess what happens for < 0. 
Let a be an arbitrary positive number. For any r > we denote 

S{r)= Yl (63) 

uiefl,0<\uj\<r 

T{r)= Yl l^r°' (64) 

ui€n,0<\u}\<r 

n{r) = #{cj G n; \uj\ < r} . (65) 

Theorem 8.13. Suppose that a = ao + a^s and that the Aharonov-Bohm vector 
potential aab is of finite type. Let the following conditions be satisfied: (a) for any 
a > 2, the sums T{r) are uniformly hounded, (6) n(r) = O(r^), (c) the sums S{r) 
are uniformly hounded for a = 2. Then, for sufficiently large 6o > 0, the Hamiltonian 
-^max(^) ^^-5 infinitely degenerate zero mode and H~^^{a.) has no zero mode. 

If for a = 2 the sums T{r) are uniformly hounded then, for any > 0, the 
Hamiltonian H^^^{a.) has an infinitely degenerate zero mode and H^^^{sl) has no zero 
mode. 

In the case when < the same claim remains true when interchanging the role 
ofH+^A^) and H-^^{a). 

Proof. Let us consider the operator H^^^{a). In view of Theorem 17. 11 we can assume 
that its zero mode, if any, has the form 

2P{z,z) = f(z)exp(^-^n^z\'^ flm''' (66) 

where Wj{z) = Wfi.{z) is the Weierstrass canonical product for the set flj (see Ap- 
pendix El) and f{z) is a nonzero entire function (cf. Lemma (8 .21) . 

Let assumptions (a), (b), (c) be satisfied. Then, according to the Borel theorem 
and to the Lindelof theorem (cf. Theorems IB. II and IB. 41 in Appendix), every function 
Wj(z) has order 2 and a finite type, i.e., |l^j(2;)| < aj exp{cj\z\'^) with some constants 
Qj, Cj > 0. It follows that for 

^>ci{i-ei) + ... + cn{i-e^) (67) 

the function is square integrable if we set f{z) = p{z) YYj=i ^ji^) where p{z) is 
an arbitrary polynomial. 

If the functions T(r) are uniformly bounded for a = 2 then < 1 (see (jlOlj) ) and, 
according to Theorem IB. 51 from Appendix, the functions Wj{z) are of minimal type 
and so the constants Cj can be chosen arbitrarily small. Consequently, the restriction 
on the field 6o > is not necessary anymore. 
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In the case of the operator H^^^{a.) we have to discuss the function 

ij{z,z) = f{z)exp(^^7T^o\z\'^ n l^^l'^"' ■ (68) 

If assumptions (a), (b), (c) are satisfied then |Wj(z)|^^^^ > aj exp{cj(9j — l)!^;^) with 
some constants aj, cj > 0. Consequently, for bo obeying (jEZj), -R > sufficiently large 
and for some c > we have the inequality > c\f{z)\'^ if \z\ > R. But then, as 

one deduces from Lemma IHUl is not square integrable, hence H~^^{a.) has no zero 
modes. 

Obviously, changing the sign at bo means that -^^^^^^(a) and H^^^{a.) interchange 
their roles in the above considerations. □ 

Remark 8.14. If conditions (a), (b), (c) from Theorem 18. 131 are satisfied then Q has 
finite density, i.e., limsup n{r)/r'^ < oo. If, in addition, the sums T(r) are uniformly 

r— >oo 

bounded for a = 2 then the density of the set Q is zero (see inequalities pOfjj) ). 

Remark 8.15. All assumptions of Theorem 18. 131 are fulfilled if every set Qj is either 
finite or a union of chains. 

Let now Qhe a lattice. Using the above introduced notation we write Qj = kj + A 
where A = {kiUi + k2UJ2] ^1,^2 G Suppose that S = lm{uiUJ2) > (S* is the 
area of an elementary cell in the lattice A). Let rjo = C,oS designate the flux of the 
homogeneous component of the field through the elementary cell of the lattice A. 

Theorem 8.16. Suppose that a = ag + a^^ and that Q is a lattice. Let b^ > 0. Then 
the Hamiltonian H^^^{sl) has an infinitely degenerate zero mode. The inequality 

n 

ilo + '^Oj < n 

i=i 

is a sufficient and necessary condition for H~^^{sl) to have a zero mode, and if it is 
fulfilled then the zero mode is infinitely degenerate. 

Let bo < 0. Then the Hamiltonian H~^^{a.) has an infinitely degenerate zero mode. 
The inequality 

n 

\Vo\ < 

is a sufficient and necessary condition for H^^^{a) to have a zero mode, and if it is 
fulfilled then the zero mode is infinitely degenerate. 

Proof. We shall start the proof from the operator H^^^. In analogy with the proof of 
Theorem 18.81 we consider the function 

n 

i^{z,z) = f{z)exp(^-^\z\')l[\a{z - (69) 
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(recall that t]q = ^qS). From Lemma f8. 71 we immediately deduce that for an arbitrary 
polynomial / it holds true that ip G L^(M^), hence ip is an infinitely degenerated zero 
mode. 

Let us now turn to the operator H^^^. According to Theorem 16.21 and Lemma f8.2[ 
the ground state of the operator H~^^, if any, has the form 

n 

where f{z) is an entire function. Using formula p26|) from Appendix [0 (where /i = 
7i/2S) we get 

n 

mzrz)\' = |/(^-)Pexp(^|.r)nexpg(0,-l)|.r) 

X f[\p{z-K,,z--^j)r^-'^ (71) 

n 

= \f{z)\' exp{c\z\') n \p{z - K„ z - 



where 

c=|(^^o+x:(%-i)j. 

The condition r^o + ^YTj=\ <'n\& equivalent to the condition c < 0. But in that case 
the membership ip G L^(M^) can be proved as in Lemma f8. 71 

Conversely, assume that c > and that there exists a non-zero entire function /(z) 
such that G L^(M^). Then from (ffT| we derive that |/(^)P < c\\'^{z^z)'\^ with some 
constant C\. Consequently, / G L^(M^) which contradicts Lemma IHTTl 

In the case 6o < one can either repeat the above reasoning or apply Corollarv l5.6l 
while noticing that {—a;} = 1 — {x} for any a; G M which is not an integer. □ 



Remark 8.17. Similarly to the case &o = (cf. Remark l8.9|l . both operators H^^J^a) 
may have localized states also when the total fiux through the elementary cell is zero. 
Suppose, for example, that 6o > 0, < < 1, < r^o + < 1 and 6^2 = 1 — r^o ~ ^i- 
Then t^o + + ^2 = 1 < 2 and the assumption of the theorem is satisfied. 

Remark 8.18. Theorem 18.161 shows that, for 6o > 0, an oscillation of the type 
"localization-delocalization" occurs after adding an Aharonov-Bohm fiux to a sys- 
tem with the Hamiltonian H~^. 
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9. Conservation of zero modes under translations 
and additions of Aharonov— Bohm solenoids. Ir- 
regular Aharonov— Bohm systems 

9.1. Translation and addition of finitely many Aharonov— Bohm 
solenoids 

Up to now we have investigated zero modes of regular Aharonov-Bohm systems (in the 
sense of the definition given at the end of Section with Theorem 18 . 1 31 representing 
the only exception. The proof of this theorem suggests that one should expect zero 
modes also in the case when the homogeneous component of the magnetic field is 
absent provided the perturbation corresponds to a (in general, irregular) "sufficiently 
scarce" system. Further we shall consider such scarce perturbations applied to systems 
of chains or lattices of Aharonov-Bohm solenoids. Before addressing this question 
we shall prove that the zero mode of the Hamiltonian corresponding to a system of 
solenoids of finite type does not disappear if a finite number of solenoids are moved 
or if one joins a finite number of solenoids to the system. 

In the foUowings theorems, a designates a potential of finite type corresponding to 
a system of Aharonov-Bohm solenoids which is determined by an array of mutually 
disjoint discrete sets, fli, . . . , fin, and by an array of fiuxes, 9i, . . . ,9n (0 < 9j < 1). Q 
designates the union Q = QiU . . .U Qn- 

Theorem 9.1. In addition to the above introduced notation let Kj = {uij, . . . ,uJn ,j} 
be a finite subset of flj and let Kj = {u[j, . . . ,uj'^^ j} be a finite subset of C such 
that the sets Q'j = {Qj \ Kj) U K'^, j = 1, . . . ,n, are mutually disjoint. If the oper- 
ator -f^max(^) ^'^^ ^ ^^^^ mode then the operator H^^^{sl ) determined by the array 
. . . , Q'^; 6i, . . . , On) has also a zero mode with the same multiplicity as that of the 
zero mode for the operator H^^^{a). 

Proof. We shall confine ourselves to the discussion of the operator H^^^{a.). The zero 
modes of this operator can be written in the form 



^{z,z) = fiz)ll\W^ 



where / is an entire function and Wj{z) = Wq. (z) is the Weierstrass canonical product 



for Qj. Then the function 



i^{z,z) = f{z)i[ [mz)r^i[ 

j=l \ k=l 



Z - UJkj 



represents a zero mode of H^^^{sl'). One has only to verify that G L^(M^). It is 
actually so because the additional singularities at the points tu^^ are square integrable 
and outside a compact set ip differs from ip hj a bounded factor. 
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This argument clearly shows that the multiplicity of the zero mode for the operator 
-^max(^') is ^ot smaller than the multiplicity of the zero mode for the operator H^^{a). 
Since the operators play an equivalent role in the assumptions the converse is also 
true. □ 

Theorem 9.2. Assume that additionally to the considered system of solenoids there 
are given a finite set = . . . ,0;^} C C not intersecting D, and a corresponding 
family of fluxes {9[, . . . ,0!^} (0 < 9j < 1). Let a! be the vector potential determined 
by the array of sets fli, . . . , Q„, D,', and by the array of fluxes 9i, . . . ,9n,0[, . . . , 9'^. 
If the operator H^^^{a.) has a zero mode then the operator H^^^^{a') also has a zero 
mode whose multiplicity is not smaller than the multiplicity of the zero mode for the 
operator H^^^^{a) . 

Proof. We shall confine ourselves to the discussion of the operator H^^{a.). The zero 
modes of this operator can be written in the form 

n 

i,{z,z) = f{z)l[\W,{z)\-'^ 

where / is an entire function and Wj (z) = Wq. (z) is the Weierstrass canonical product 
for It turns out that the function 

n m 

i^{z,z) = f{z)i[ mz)r^ i[\z - . 
j=i k=i 

is a zero mode of H^^{a'). One has only to verify that ip e L^(R^). This again follows 
from the fact that the singularities at the points are square integrable and the 
function ip differs from by a bounded factor outside a compact set. □ 

9.2. Additional notation 

Up to the end of the current section, flj (j — 1, . . . ,n) designates either an array of 
mutually disjoint chains or an array of mutually disjoint lattices, ftj — Kj + Aj where 

Aj is a Bravais lattice of rank 1 or 2 and Kj = {nij, . . . , K-m^j} is a finite set. To each set 
Qj we relate an array of Aharonov-Bohm fiuxes Qj = {9kj)i<k<m ■ By Vt wc denote the 
union 1] = fiiU. . .Ufi^. In addition, we shall consider another array of discrete subsets 
in the plane, 0'^,...,^^, whose members are mutually disjoint as well as disjoint 
with the sets Qi, . . . , and we relate to these additional sets an array of fluxes 
®' — {^j)i<j<m- Finally, we consider a discrete set C whose points are supposed to 
be removed from Q. Set Q' = Q'qUQ'iU . . -Ufi^. Furthermore, Wj{z) is the Weierstrass 
canonical product for the set fl'j, Vkj{z) is the Weierstrass canonical product for the 
set VIq n {nkj + Aj). By r' we denote the convergence exponent of the set Q,' and by p' 
its genus. The symbol n'{r) designates the number of points of the set Q' contained 
in the disc \z\ < r. The symbol a' designates the vector potential for the perturbed 
system of solenoids determined by the discrete sets f2i \ JIq, . . . , 0„ \ JIq, il'^^ U . . . U fl'^, 
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and by the array of fluxes obtained by concatenating the arrays Gi, . . . , and 0'. 
Let us note that we still assume that < 6*^ < 1, < 6'^- < 1. 

Thus the set fl' representing the total perturbation of the original set Q need not 
be finite nor regular. On the other hand, we shall always suppose that the set Q' is 
sufficiently scarce by imposing restrictive assumptions on its genus p'. 

9.3. Addition of solenoids to a union of Aharonov— Bohm chains 

In Theorem 18.41 it has been shown, roughly, that if is a finite union of chains then 
the Hamiltonians H^^^[Q) have infinitely many zero modes. Below we show that this 
property survives provided the perturbation Q' is sufficiently scarce and at least too 
chains are not parallel. 

Theorem 9.3. Let . . . , fi„ be chains whose union is a uniformly discrete set. Sup- 
pose that among the chains there are at least two which are not parallel. Furthermore, 
suppose that the genus of fl' fulfills p' = 0. Then the Hamiltonians H^^^{sl') have 
infinitely degenerate zero modes. 

Proof. In virtue of Lemma fA. II one can assume that every chain Qj = Kj + Aj, with 
Aj = {kuj] k G Z}, ujj 7^ 0, is contained in a line Lj and that different chains are 
contained in different lines. We shall suppose, without loss of generality, that Li and 
L2 are not parallel and that Li coincides with the real line M, with G i^i. 
Let us consider a function ip of the form 

n m 

i,{z,z)=f{z)\[g,{zrz)X{\W,{z)\'-''^ (72) 
i=i k=i 

where 

9j{z,z) = Yl 

k=l 

To prove the theorem it suffices to find an infinite, linearly independent family of 
entire function f{z) such that ip G L^(]R^). 
Let us show that the functions 

fiz) = (74) 

with sufficiently small a > suit this condition. To this end we shall need the following 
lemma. 

Lemma 9.4. Assume that ai, 0:2, cti, 02 G C fulfill axa^ 7^ 0, axja^ ^ M. Then for 
every e > there exist constants c, Ci, C2, 71, 72 > such that the inequality 

cie^il^l < I sin(ai(;z - oi)) sin(a2(2 - 02)) I < C2e^2l^l (75) 

holds true whenever \z\ > c and the distance from z to the lines Li = a^^M. + ai and 
L2 = a^^M + a2 is greater than e. 



sm 



nlz 



Kkj) 



kj\ 



(73) 
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Proof of Lemma \9.4\ Let L be a line written in the form L = a ^M+a where a, a G C, 
a ^ 0. Then inequahty (jKn|) imphes that 

^^^^^ < |sin(a(z-a))| < el"!'^ (76) 

where d is the distance from the point z to the hne L. Actually, set a = \a\e^'^. Then 
I ImQ;(2; — a)\ = \a\d where d is the distance from e*'^(z — a) to M. But, at the same 
time, d is the distance from z to e~'^'^'R + a = a~^M. + a. From ()76|) we deduce that for 
every e > one can find a constant c > such that 



sm 



(a(z-a))| > cel"!'^ (77) 



whenever the distance d from 2 to L is greater than e. Actually, it suffices to choose 
c= (l-e-l"l^)/2. 

Let us now reconsider the lines Li and L2 from the lemma. Since ai/a2 ^ M the 
lines intersect each other in a point f G C. Let us denote by ip the angle between the 
lines Li and L2 {0 < ip < n), and by 9 the angle between the vector z — v and the 
line Li. Then the distances di and d2 from the point z to the lines Li and L2 are 
respectively equal 

di = \z — v\ \ sin(6')| , d2 = \z — v\ \ sin {9 — ip) \ . 

From dZni) we get 

I sm{ai{z - ai)) sin(a2(2 - ^2)) | < el^^l^^+'^^l'^^ < e'^^^d^il'I^^DC^i+'ia) , 
Notice that di + d2 < 2\z — f | < 2\z\ + 2\v\, hence 

I sin(Q;i(2; — ai)) sin(a;2(-2 — 02))! < C2e'''^'^ 

where 

72 = 2max(|ai|, \a2\) , C2 = e'"''^^ , 

(this is true for any z G C). 

Using (f77j) we can relate to every e > a constant c > such that if di,d2 > s 
then 

I sin(a(z - ai)) sin(a(z - 02)) | > c2e"^''^(l"il'l"2l)("'i+'^2) . 
On the other hand, 

di + d2 > 1-2 — f |(sin^ 6* + sin^(^^ — (y?)) = |z — t>|(l — cos(/9cos((y9 — 26*)) 
> \z — v\{l — I cosy^l) > \z\{l — I cosy^l) — |t>|(l — | cosy?!) . 

From here we deduce that the inequality 

I sm(^a{z — oi)) sin(a(2; — 02)) | > Cic'''^'^' 

holds true for 

ci = c2e-"^i'^(l°il'l°2l)l'''l(i-lcosv.|) ^ _ min(|ai|, |a2|)(l - | cos^\) . 
This proves the lemma. □ 
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Proof of Theorem \9.'^ (continued). Let us return to the proof of Theorem 19.31 From 
the assumptions of the theorem [p' = 0) it follows that the functions Vkj{z) and Wj{z) 
are of growth (1, 0) (see Appendix Ej) • This implies that for every e > there exists a 
constant > such that 



n rrij 



j=i k=i 



< Ceexp(e|2:|) . 



(7J 



1=1 



Let Pj be a strip with border lines parallel to Lj and containing Lj in its interior, 
and let Q be a sufficiently large disk centered at 0. In virtue of Lemma I9.4| the disk 
Q can be chosen so that, for 2; ^ Q U Pi U P2, the inequality 



{zrz)\<c^\f{z)\\~gf{zrz)\\~9i{zrz)\\{\~9, 

i=3 



7.A, 



Z,Z) 



(79) 



holds true with some constant ci > 0. Here we have set 

k=l ^ 

and 



~9fiz,z) 



sm 



( 



niz - Ki,- 



n 

k=2 



sm 



with < (3j < 6ij. 

On the other hand, one can choose Q so that, for z ^ P2\Q,we have the inequality 
\z\ < c'di where di is the distance from z to the line Li, and c' > does not depend 
on z. Then from ()77j) we deduce that Q can be replaced by a larger disk such that the 
inequality 



z,z] 



0) 



\i;{z,z)\<c\\f{z)\\~g^{zrz)\\{\~9', 

with a constant c'^^ > 0, holds true for z G P2 \ <5- An analogous assertion is true when 
interchanging the strips Pi and P2- 

Finally, for any choice of the disk the inequality 



\i;{zrz)\<c'i\f{z)\\{\i 



'i{z,z) 



^1) 



holds true in the interior of Q. Formulas (|7n|) . (jHOj) and (jHT|) make it possible to 
complete the proof by arguing in the same way as in the proof of Theorem 18 .41 □ 

In the case when all chains are parallel we have a somewhat weaker result. 

Theorem 9.5. Let fii, . . . , f2„ he parallel chains whose union is a uniformly discrete 
set. Assume that the convergence exponent of Q' satisfies either t' < 1/2 or r' < 1/2 
and n'{r) = o(r^/^). Then the Hamiltonians H^^^^(sl') have infinitely degenerate zero 
modes. 
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Remark 9.6. Under the assumptions of Theorem 19.51 it holds true that p' = but 
this equahty does not imply the assumptions of the theorem. 

Proof. In virtue of Lemma lA.ll one can assume that every chain flj = Kj + Aj, 
Aj = {kuj] /c G Z}, is contained in a line Lj, with different chains being contained in 
different lines, and that all lines are parallel to the real axis. Hence one can assume 
that ujj > for all j and that all lines Lj are contained in a half-plane Im z > a where 
a > 0. 

Let us consider a function i/) of the form 

n m 

^(z, -z) = f{z) n 9,{^, -z) n \Wk{z)\'-'''' (82) 

j=l k=l 



where 



9j{z,z) = Yl 



k=l 



sm[ 



\VkAz)r^- (83) 



To prove the theorem it suffices to find an infinite, linearly independent family of entire 
functions f{z) such that G L^(]R^). Let us show that in this case the functions 

fiz) = . , , — ; ^=v , (84) 

sm(^7razj ^m.[^\/ —t{ az) 

with sufficiently small a > 0, will do. Here the function /(z) is well defined and 
analytic in the upper half-plane provided the usual branch of the square root has been 
chosen. If Imz > then ^/z^J—z = —iz. Since sm{^/z)/^/z = 1 — z/3\ + ... is in fact 
an entire function, also the function f{z) extends as an entire function. 

We start the verification from several preliminary observations. The first one follows 
from inequality ()50|) . 

(A) For every e > there exists c > such that 

I sin(v^)| > - expd^;!"*"/^) for > c, £ <\ aigz\ < vr, 
3 

I sin(V— -z) I > ~ Gxpd^;!"*"/^) for > c, < | arg2;| < vr — £. 
3 



Suppose that n G N and < 5 < 7r/4. Let us denote 

Bn{5) = {2; G C; \\fKZ — -nn\ < 5} . 
(B) For n 7^ m the sets -B„(5) and Bm{6) are disjoint. 

Actually, suppose that n > m. If 2 G Bn{S) fl Bm{S) then there exist u,v G C, 
\u\, \v\ < 6, such that nz = {im + -u)^ = (vrm + f)^. Hence 7r^(n — m){n + m) = 
271 {mv —nu) +v'^ —u"^ . At the same time it holds true that 7r^(n— m)(?2+m) > TT^{n+m), 
\27r{mv — nu) + f ^ — < 27r6{m + n) + 25^ < 47r5(m + n) < n'^in + m). 
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Set 

Q{e) = {z eC; \ sin(v^)| < e}. 

Let us denote by Qni^) the connected component of the set Q^e) containing the 
point nn^, and by U{e) the connected component of the set {w E C; \ siia{w)\ < e} 
containing 0. Observe that (B) imphes the following claim. 

(C) For sufficiently smalls > the connected components Qni^) are mutually disjoint. 
To complete the proof we shall need the following two lemmas. 

Lemma 9.7. For every S > there exists a constant c > such that 

< cmax(|v^7r^| + 1, | sin(2)|) 



sm(z) 



sm 



sin^(x) + sh^(?/), with z = x + iy, 



on the strip \lmz\ < 6. 

Proof of Lemma \9. 7| The equality |sin(2;)p 
a;, 1/ G M, implies that | sin z\ < (x^ + cjy'^Y^'^ < Cs\z\ for | Im2;| < S where we have 
set cs = sh((5)/5. 

Let us choose e > small enough so that the sets Qni^) are mutually disjoint, 
\^/^Tz\ > nn — 1 for z E Qni^) and | sin(w)| > \w\/2 on U{e). For z G C \ Q{e), the 
desired inequality is valid with c = e~^. If 2; G Qni^) then a/ttz — tttt, G U{e) and 
therefore 



I sm[y/7rz)\ = 
Furthermore, if | Imz| < 6 then 

I sinfz)! = 



sin (y/nz — nn)\ > \y/nz — 7m\/2 . 



sm(z — Tin )| < cs\z — im 



Consequently we have, for z G Qni^) and | Im2;| < 5 

t2 



smu 



<2C5 



7rn 



vrn 



TT 



smi^yvrz 
This proves the lemma. 

Lemma 9.8. For any b > there exists e > such that 



< — (Iv^l + vrn) < — (Iv^l + 1) . 

~ TT 



□ 



smu 



sm 



exp(^ — b\z\^^'^) dxdy < 00 



Qie) 



where one chooses the principal branch of the square root on C\ M_. 

Proof of Lemma \9. ^ We choose e small enough so that Claim (C) is true. In the 
integral 

exp( — dxdy 





sin(2;) 




sin(A/7ri) 



(3n(e) 
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TTZ — 7rn, i.e., z = {w + 7rn)^/7r, where w = u + iv, 



we apply the substitution w = 
u,v eR. Since dx A dy = (4/7r^)|w + nn\'^du A dv and Qni^) is mapped onto U{e) we 
get 



C/(e) 



sin ( 2nw + — 



sm Uf 



|w + vml exp —\w + 7m\\dudv. 



One can assume that e is sufficiently small so that | sin(iy)| > 11^1/2 on U{6). Further- 
more, there clearly exists a constant c > 0, depending on e but independent of n, such 
that I sm{2nw + {w'^/n))] < csh(2?T,|w|) on U{e). Thus we get 



16c2 



sh^(2n 


\w\) 




w\ 


2 



|w + 7™| exp 



\w + TTn\ I dudv 



TT 



U(e) 



By modifying the constant in front of the integral we can simplify this inequality. 



In < c'n^ exp ( — & n) 



sh^(2n 


\w\) 




w\ 


2 



dudv . 



U{e) 



Here again the constant c' > does not depend on n. If sup{|w|; w G f/(£)} < by/n/i 
then 



n=0 



5^4 <c' // J^n^exp 



sh\2n\w\) 



dudv < oo - 



U{e) 



n=0 



\W\ 



This proves the lemma. 



□ 



Proof of Theorem \!J.,^ (continued). Let us return to the proof of Theorem 19.51 We 
denote Ai = {z E C; Re z > 0}, A2 = {z E C; Re z < 0}, and we shall show that 
ijj e L'^{Aj), j = 1, 2. More precisely, we shall prove only the membership ip G L^(Ai), 
the property G L'^{A2) can be shown analogously. We split the set Ai into a union 
Ai = PiU P2 with P^ = {z e Ai; < Im^ < 6}, P2 = Ai\ Pi. The bound b is chosen 
so that the strip Pi contains the set Q in its interior. 

From the assumptions it follows that the functions Vkj{z) and Wj{z) are of growth 
(1/2,0), i.e., for every e > there exists a constant > such that 



maXj^k{\Vkj{z)\, \Wk{z)\} < Csexp{e\ 

1 



(cf. Appendix IB]) . Let us denote 



G{z,z) 



Hence 



where 



i){z,z) 



j=i k=i 



^5) 



smiaz) 



sin{y^7raz 

rrij 



G{z,z)Y[9j{z,. 
i=i 



k=l 



sm 



7r[Z - Kki 
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According to Claim (A) there exist constants Ci,C2 > such that 

\G{z,z)\ <ciexp{-C2\z\^/^), (86) 

for all z ^ Ax. 

With the aid of Lemma 19.71 and formula (jHBj) we can estimate the function t/^ on 
the strip Pi, 

n 

\i){z,z)\ < c'dVvra^l + 1) exp(-C2|2;|^/^) JJ ^^(2;, ^) . 

The singularities of ^/^ in Pi are square integrable and therefore, similarly as in the 
proof of Theorem 18.41 we obtain 

(l)^GL2(Pi). 

Since C Pi, on P2 fl a~^Q{e) we have the estimate 



[zrz\\<d' 



sm 0:2; 



exp(-C2|z 



1/2N 



sin(-y7rai) 

If e is small enough then Lemma [9.81 implies that 
(2) i) e L\P2na-^Q{e)). 

Finally, for 2; G P2 \ a^^Q{e) we have | sin(-^7ra;2;)| > e and hence the inequality 

n 

\iIj{z,z)\ < c"'| sin(az)| 

i=i 

holds true. We conclude, similarly as in the proof of Theorem 18. 4| that if 

n rrij 

j=i k=i ^ 

then 

{3)^eL\P2\Q^{e)). 

This concludes the proof of the theorem. □ 

9.4. Addition of solenoids to an Aharonov— Bohm lattice 

The case when the Aharonov-Bohm fluxes are arranged in a lattice Q has been dis- 
cussed in Theorem 18.81 It turns out that that for a scarce perturbation the result 
stated in the theorem is still true. 

Theorem 9.9. Suppose that Q is a lattice, i.e., Qj = kj + A where A is a Bravais 
lattice of rank 2. Suppose further that the genus of the set Q' fulfills p' < I. Then the 
Hamiltonians H^^^la') have infinitely degenerate zero modes. 
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Proof. Let W{z) be the Weierstrass canonical product for the set Q', let g' be the 
growth order of W{z), and let r' be the convergence exponent of Q'. Then g' = t' < 
p' + l<2 (see Appendix E)). If ^' < 2 then W{z) is of growth (2, 0). If ^' = 2 = + 1 
then, by Theorem IB. 51 (h). W{z) is of minimal type. Consequently, also in the latter 
case W{z) is of growth (2,0). This means that for any c > there exist a > and 
R> such that for all z, \z\ > R, it holds true that 

\W{z)\ <a exp{c\z\'^). (87) 

The same observation is clearly true for any subset of Q'. In particular, the functions 
Vj{z) and Wi{z) are of growth (2,0) and obey estimates similar to (|S7|). 
Zero modes of H^^^{a.') are gauge equivalent to functions of the form 

^(z, -z) = f{z) n - ^^r'^ n n ■ 
j=i j=i j=i 

From Lemma 18.71 and from the estimate ()87|) with sufficiently small c > it follows 
that ip is square integrable if f{z) is an arbitrary polynomial. □ 



9.5. Irregular systems of Aharonov— Bohm solenoids 

All the preceding results were concerned with Aharonov-Bohm systems Q with bounded 
density, i.e., for which limsup^^g^ n{r)/r'^ < oo (cf. ( Here we show that zero 
modes may occur also in systems with infinite density, more precisely, in systems for 
which liminf n(r)/r^ = oo. Moreover, we shall present examples of systems Q with 

r— >oo 

arbitrarily large convergence exponent tq. Let us fix a natural number N > 2 and set 

= {e^^^^/^ m^/^; m G N , = 0, 1, . . . , 2iV - 1} . 

Obviously, the convergence exponent r for the set Qn equals A^. In particular, for 
A^ > 2 we have lim^^oo n{r)/r'^ = oo. Let 6 be an arbitrary number from the interval 
]0, 1[. Then the vector potential a of the Aharonov-Bohm system determined by the 
couple {Qiy,9) reads sl{z,z) = 6 sgrad ln(|iy(2;)|) where 



W{z) 



sin(7r2;^) 



Theorem 9.10. The Hamiltonians H^^J^a) have infinitely degenerate zero modes. 
Proof. It is sufficient to show that for < a < 7r6' the function 



sin(az^) 



is square integrable. Set 
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Then 

p{z, z)]"^ dxdy = N // \'ip{z, z)]"^ dxdy 



s 



and therefore it suffices to verify that 

Ij{z,z)\'^ dxdy < oo . (88) 



s 

Let us make a substitution of the integration variable in (jHH|l . w = z^ where w = u+iv, 
M, f G M. Since du A dv = N'^\z\'^'^~'^dx A dy we can rewrite the integral as 

jj(z, z)]"^ dxdy = JJ j — ^\ sm{aw)\'^\ sm{Tcw)\~'^^ dudv (89) 
where 

p = 4-26-2^--^. 

N 

Since 2 — 29 — j3 > —2 and /3 > 1 one can show that the integral ()89p is finite using a 
reasoning as in the proof of Theorem 18.41 □ 



Appendix A. Lattices 

Here we collect basic definitions and some auxiliary results about lattices. Let E be 
a finite-dimensional real Euclidean space with dimension d. A discrete subgroup A of 
the additive group E is called a Bravais lattice. For any Bravais lattice A there exist 
linearly independent vectors cji , . . . , G -E such that 

A = ZtUi + Zco'2 + . . . + IjUJr ■ 

The array {ojj)i<j<r is called a basis of the Bravais lattice A. The integer r does not 
depend on the choice of basis and is called the rank of the lattice A. To every basis 
{i-Oj)i<j<r one relates the elementary cell F, F C E, formed by all x E E whose 
orthogonal projection x' onto the linear span L of the lattice A has a decomposition 

x' = tiLUi + tiU2 . . . + trUJr , 

with <tj < 1 for all j. li r = d (the dimension of the lattice A is maximal possible) 
then F is a convex parallelepiped. In the opposite case F is even not bounded. 

A non-empty discrete subset Q G E is called a crystal with the Bravais lattice A if 
it is invariant with respect to the action of A on and has a finite number of orbits. 
Obviously, every crystal fl can be written in the form Q = K + A where K G E is 
a finite set whose number of elements equals the number of orbits. Without loss of 
generality we may assume that K G F. Conversely, every set of the form Q = K + A 
is a crystal. If = 1 then the crystal is called mono- atomic or simple. In the general 
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case when \K\ = n the crystal Q is called n-atomic. If r = 1 then Q is called a chain (a 
simple chain if in addition \K\ = 1). If r = d then Q is called a lattice (more precisely, 
a crystal lattice) in the space E. In other words, a crystal is such a discrete subset 
Q G E whose group of parallel translations acts co-compactly on E. 

Let us note that in our definition we do not exclude the case r = 0. If so then 
A = {0} and a crystal with the Bravais lattice A is simply a finite subset of E. 

It is worth of noticing that a crystal Q is always a uniformly discrete subset of 
E. This means that there exists a constant c > such that — uj"\ > c whenever 
u', uj" en,iu'y^ Uj". 

We shall need the following lemma. 

Lemma A.l. Assume that dimE = 1 and that fii , . . . , f2„ are chains in E. The 
union Q = Qi U . . . U Qn is a chain if and only if Q is a uniformly discrete set. 

Proof. We only need to prove that this condition is sufficient. Moreover, it suffices to 
consider the case n = 2. The general case then follows by mathematical induction. Let 
us write 

nj=Kj+Aj (j = l,2), 

where Aj is the Bravais lattice of Qj. Let us identify E with M. Then Aj = TjUj, with 
ujj > 0. 

We shall show that the number p = uj^juj^ is rational. Actually, in the opposite 
case the set Zcui + Za;2 would be dense in M. Let us choose k\ G K\ and 1^2 & K2. We 
can find a sequence nkUJ2 — rrikUJi {rrik ,nk E Z) converging to Ki — K2 and such that 
/^i — /«2 7^ nkUJ2 — fnkUJi for all k. Obviously, this contradicts the assumption that the 
set Q is uniformly discrete. 

Hence p = N/M, with iV, M G N. Then Mui = Nuj2 and therefore VL is invariant 
with respect to the lattice A with the basis vector Muji. It is easy to see that the 
number of orbits of the group A in f2 is finite. But this means that is a chain. □ 

Remark A. 2. For dimii^ > 1 an analog of Lemma fA. II is false. Actually, already for 
dim E = 2 it can happen that a union of two simple lattices is uniformly discrete but 
not a lattice. This is demonstrated by the following example. Let E = M^. Let Ai be 
a Bravais lattice with the basis ui = ei, uj2 = 62, and let A2 be a Bravais lattice with 
the basis uj[ = V2ei, ^'2 = ^2- Consider the crystal lattices fii = Ai and ^2 = + ^2 
where k = (l/2)e2. Obviously, f2 = f^i U 1^2 is a uniformly discrete set. Let A be a 
group of parallel translations acting on Q. Suppose that the number of orbits of the 
group A in n is finite. Then there exist n , m G Z, n 7^ m, such that nui and mui 
belong to the same orbit. Hence kui G A for some /c G Z, A; 7^ 0. But k + kui ^ Q. 

Appendix B. Auxiliary results from the theory of 
analytic functions 

Here we recall some results from the theory of analytic functions that are necessary 
for our presentation, for the details see ll()2| llOfjj . For an entire function / we set 
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Mf{r) = M{r) = max \f{z) \ = max |/(^)| . (90) 

\z\=r \z\<r 

The order (more precisely, the growth order) of an entire function / is the number 

Qf = Q^ inffa; 3Ra >Oyr>Ra M{r) < exp(r")} , (91) 

or, equivalent ly, 

y ln(ln(M(r))) 
^ ^r-.r ln(r) ■ ^'2) 

Let us note that if f{z) = Y^'^=o ^^^^ th.%TL 

\w\Tl] 

If ^/ < oo then one says that / is a function of finite order. For a function of finite 
order q the number 

^/ = ^ = m.i{K > 0; ^Rk >Q\/r> Rk M{r) < exp{Kr^)} (94) 

is well defined and it is called the type of the function f. The type can be equivalently 
defined by the formula 

r HMjr)) 
qf — hm sup — ^ . (95) 

r— >oo r^ 

Moreover, qf can be expressed in terms of the Taylor coefficients a„, 

i^fegf)^/^ = limsup n^^^ |a„|^/" . (96) 

n— +00 

If ^/ < OO and = then / is called a function of minimal type. 

A function / of finite order g and of finite type q obeys the estimate 

\f{z)\<eMi^ + ^)m (97) 
for arbitrary £ > and \z\ greater than a constant R^ > 0. Conversely, if the estimate 

\f{z)\<cexp{q,\zn (98) 

is fulfilled with a constant c > then the function / has both a finite order and a finite 
type, and it holds Qf < gi and qj < qi. A couple of numbers (gi.qi), < gi,qi < oo, 
determines the growth of a function f{z) if gf < gi and gj = gi implies qf < qi. 
Functions of growth (l,^i), with < oo, are said to have exponential growth. 

The order and the type of an entire function on one side and the distribution of 
its zeroes on the other side are deeply related. Let us arrange all nonzero elements of 
a discrete set O C C in a sequence = {u!k)k>i which is ascending in the absolute 
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value and ascending in the argument (0 < arg z < 2tt) in the case of equal absolute 
values. The convergence exponent of the set Q (or of the sequence Q^,) is the number 



Tn = T = ini ■{ a > 0; } ^ — ^ < oo ^ , (99) 

or, equivalently. 



ln(A;) 

tq = limsup — — - . (100) 



If tq is finite then the number 



oo 



max < n G N; V = oo> if is infinite, 
Pn^P={ I J (101) 

—oo if n is finite, 

is well defined and it is called the genus of the set Q (or of the sequence fl^). For 
tq = cxD we set pq = oo. 
For r > we set 

nn{r) = n{r) = #{c<j G Q; \uj\ < r} . (102) 
For a non-empty set Q the formula 

ln(n(r)) 

Tn = limsup , , , , (103) 
r^oo ln(r) 

shows that that the convergence exponent of an discrete set characterizes its density 
[TOHI Theorem 2.5.8]. 

On the other hand, let f2/ = be the zero set of an entire function /. Then the 
following fundamental inequality is valid (the Hadamard theorem, see for example 
jTM Theorem 2.5.18]): 

Tn<Qf- (104) 
If Qf is not an integer then tq = Qf |106t Theorem 2.9.1]. From ()104|) we deduce that 

n(r) = 0(r^+^) (105) 

for any e > 0. If > and <jj < oo then a stronger estimate is valid |106t Theorem 
2.5.13]: 

L = lim sup r^^n{r) < CQf (^f , 

r^oo 

I = lim inf r~^n{r) < Qf^f- (106) 

Moreover, if > then Le^/^ < Qji^f, in particular, L + I < Qf'^j- 

For Qf < oo it can never happen that nn{r) = o(r^~^) jl06t Theorem 2.9.3]. The 
following theorem is due to Lindelof (see jl06| Theorem 2.9.5 and Theorem 2.10.1]). 
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Theorem B.l. (1) Assume that q = Qf < oo is not an integer. An entire function 
f{z) is of finite type if and only if nuir) = 0{r^), and it is of minimal type if and 
only if n^{r) = o(r^). 

(2) Assume that qj is a positive integer. The function f{z) is of finite type if and only 
if nn{r) = O(r^) and the sums 



S{r)= 

0<K|<r 



are bounded. 



An entire function with simple zeroes is determined by its zero set up to an multi- 
plier e^^^^ where g{z) is an entire function. Furthermore, for an arbitrary discrete set 
f2 C C there exists an entire function f{z) with simple zeroes whose zero set coincides 
with Q. Denote by E(u,p) the Weierstrass canonical multiplier, with m G C and with 

pen, 



E(u, p) = (1 - u) exp M + — + ... H 

V 2 p 

(by definition, E{u, 0) = 1 — u). Let VL^ = {u!k)k>i be, as above, the sequence formed by 
all nonzero elements of the set Q appropriately enumerated. Let us denote hj xn = X 
an integer which is equal 1 if G and in the opposite case. The Weierstrass 
canonical product associated to Q is by definition an entire function Wn{z) defined by 
the infinite product 

oo 



l[E{z/uJk,Pu) (107) 

k=l 

if the convergence exponent of Q is finite, and by the infinite product 

oo 

l[E{z/uJk,k) (108) 



z-^ 

k=l 



in the opposite case. 

Theorem B.2 (Weierstrass, Hadamard). The infinite product defining the func- 
tion Wn{z) converges absolutely and locally uniformly. Consequently, W^{z) is an 
entire function and its zero set coincides with Q. Moreover, the zero set of an entire 
function f{z) with simple zeroes only equals Q if and only if the function f{z) is of 
the form 

fiz) = e^(^) Wniz) (109) 

where g{z) is an entire function. The growth order of the function Wn{z) equals the 
convergence exponent of the set Q. 

Moreover, the following theorem is true. 

Theorem B.3 (Hadamard). If the function f from relation MU9ji has a finite order 
Qf then g{z) is a polynomial with a degree not exceeding [gf]. 
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Theorem B.4 (Borel Theorem). Conversely, if tq < oo and f{z) is a function 
written in the form MO 9^) where g{z) is a polynomial of degree n then f{z) has a finite 
order Qf = max(r, n). // either tq < n or the series YlT=i l^fcl""" ^■^ convergent then 
the function f{z) is of finite type. 

The genus of a function f{z) having the form ()109|) . where g{z) is a polynomial 
of degree n, is the integer qf = q = max(n,pn). The following theorem is a useful 
completion of Theorem IB. II due to Lindelof |106| Theorem 2.10.3]. 

Theorem B.5. Under the assumptions of Theorem \B. 1\ let Qf he a positive integer. 
A function f{z) written in the form MOf^) . where g{z) is a polynomial, is of minimal 
type if and only if one of the following conditions is satisfied: 

(a) n^{r) = o{r^), pn = Qf, and 

oo 

^uj~^ = -Qfao 

k=l 

where is the coefficient (possibly vanishing) standing at z^ in the polynomial g{z), 

(b) Pn = Qf — ^ and ao = 0. 

In particular, if qf < Qf then f{z) is of minimal type. 

We shall also need the following particular case of the Mittag-Leffler theorem (see 
II.7.3.2] or [TU?). 

Theorem B.6. For an arbitrary discrete subset Q of the complex plane C and for an 
arbitrary sequence of complex numbers {6uj)ujen there exists a meromorphic function 
M{z) obeying the following conditions: 

(1) M{z) has only simple poles, 

(2) the set of poles of the function M{z) coincides with Q, 

(3) the residuum of M{z) at the point oj equals 9^. 



Appendix C. The Weierstrass cr- function and related 
functions 

In our approach an important role is played by the order q and by the type ^ of the 
Weierstrass cx-function (t{z), 



a[z) = a[z;uji,uj2 

ioeA\{o} 

It is easy to see that the convergence exponent of any lattice in the plane equals 2. 
Actually, the series 

+0O 



E' 1 r (110) 



ni,n2=— oo 
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(the dash indicates, as usual, that the summand with indices ni = n2 = is omitted) 
converges if and only if a > 2. Hence, by the Borel theorem, q = 2. Since for a = 2 the 
series diverges the Borel theorem does not say anything about the type of the 

function cr(z) (apart of the fact that it is finite). The type of this function has been 
found in the general case by A. M. Perelomov [^3]. In order to make our presentation 
self-contained we reproduce below some details from his derivation. 
Let us start from recalling the notation 

^. = 2C(f) (111) 
and the fact that the cr-function is quasi-periodic in the following sense: 

a{z + uj,) = -a{z)e^^[iqj[z+'^y^ . (112) 

Recall also that S = Im(ti;iCo'2) designates the area of the elementary cell. 
Lemma C.l ([95j). The function \(7{z)\'^ can be expressed in the form 

|cr(z)|^ = exp{h'z'^ + uz^ + 2iizz)p{z, z) (113) 
where p is a A-periodic function, 

^ = ^(^I'^S - r/2(^l) ' = ^ • (114) 

Proof. From ()112|) we obtain 

|a(2; + a;,)|2 = |a(2;)|2exp(2Re(r/,(z + ^))) . (115) 

On the other side, the function p defined by equality ()113|) . with i/ G C and /i G M, is 
periodic if and only if it holds 

\a{z + ujj)\'^ = exp(2uzujj + 2uzujj + 2p{ziUj + zcuj) + i/cul + ucu'^ + 2pujjUJj) \<y{z)\^ . 

(116) 

Comparing ()115p to ()116|) and taking into account the equality 



2Rer],(^z + ^) = r]jz + r],z + + ^ 
we arrive at the system 



uuj + puj = -r]j , i = l,2, 

uuj'^ + z/cD| + 2pujjU!j = ^{r]jUJj + VjUj) , j = 1, 2 . 



(117) 



The first couple of equations in ()117|) gives 



^ = 7^— = > /^=o— = • (118) 

2 UJ1UJ2 — UJ1UJ2 2 UJ1U2 — UJ1LU2 
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Since uj\uj2 — ^\'^2 = — 2iS' and in virtue of the Lagrange identity 

?7icj2 - ri2UJ\ = 2m (119) 

we find that relations (jllSj) and ()114|) are equivalent. Using ()118p and the fact that /i is 
real one can check that the second couple of equations in ()117|) is satisfied identically. 

□ 

Lemma C.2 ([95]). The type of the function cr(z; tui, ^2) is given by the equality 

= \u\ + fi = ^ {\r]iLJ2 - r]2i^i\ + 271) . (120) 

Proof. Let us rewrite ()113|) as follows, 

|(t(z)P = exp((z/ + z/)(x^ - y^) + 2i{u - v)xy + 2/i(x^ + y"^)) p{z, z) . (121) 

The quadratic form occurring in the exponent, 2Re(i/)(x^ — y^) — 4Im(z/)xy, can be 
diagonalized with the aid of a rotation of the coordinate system. Eigenvalues of the 
corresponding symmetric matrix are Ai = — A2 = 2|z/|. Set e = e^'^ where is the angle 
of the rotation. Since the quadratic form + is rotationally invariant we have 

\a{ez)\ = exp {{ft + \u\)x^ + {ft - \u\)y^) p^/\ez,ez). (122) 

Owing to the periodicity of the function p it holds true that 

max \<j{z)\ = max |cr(ez)| < cexp ((/i + |i^|)r^) . (123) 

|2|=r |2|=r 

Consequently, <j < |z/| + /i. 

To show the opposite inequality it suffices to construct a sequence z^ such that 
\zk\ 00 and 

\cr{ezk)\ > cexp((/i + \u\ - 6k)\zk\'^) , (124) 

where Sk I and c > is a fixed constant. First we note that, by the uniqueness 
theorem for analytic functions in a real variable, there exists a point zq such that 
p{zo, zq) 7^ 0. Then there exists c > such that \p{z, z)\ > c on a neighborhood V of zq. 
This gives the choice of c. Further we consider the canonical mapping h : — > M^/A. 
Two cases are possible: either the image h{zo + M) is a closed curve in the torus 
T = R^/A or this image is dense in T. In the former case there exists a sequence 
Afc e M such that Xk —>■ 00 and 

hizo + Xk) = h{zo) , (125) 

in the latter case condition ()125|) should be replaced by h{zo + Xk) h{zo). In the 
both cases condition ()124|) holds true with Zk = zq + Xk- □ 

Following jnS] we introduce the function 

a{z) = e-"^V(z). 

Lemma fC. II implies the equality 

\a{z)\'^ = exp{2p\z\'^)p{z,z). (126) 
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Lemma C.3 ([95]). Let f{z) be an entire function whose zero set coincides with 
A = IjjJi + 'Luj2, with all zeroes being simple. Then the order Qf is at least 1, gj > 2, 
and if Qf = 2 then the type ^/ is at least /i, ^/ > yU = 11/28 . Moreover, in the case of 
the function a{z) the minimal values are achieved both for the order g and the type 
i.e., Qa = 2 and q^. = fi = tt/2S. 

Proof. Since the function a{z) is expressed as a Weierstrass canonical product its 
order equals the convergence exponent ta = 2. Let us consider the entire function 
f(z) = e-"^V(;z), with a G C. Then 

= exp(^2Re(z/ - a){x^ - y^) - 4Im(z/ - a)xy + 2fi{x^ + y^)^ p{z,z) . (127) 

It is clear that the order of the function f{z) equals 2, and similarly as in the proof of 
Lemma Rl2| the type of / equals |z/ — a| + /x. Obviously, the smallest type (namely, /i) 
is achieved for a = u. In particular, the function cr(z) is of order 2 and its type equals 

Conversely, suppose that the zero set of an entire function f{z) coincides with A 
and that all zeroes of f{z) are simple. Since t\ = 2 the Hadamard theorem implies 
that Qf > 2. Suppose that gj = 2. We can write f{z) in the form f(z) = e^^^^a{z). By 
the Hadamard theorem, g{z) = az"^ + bz + c. If a = then the type of f{z) equals the 
type of cr{z), if a 7^ then the type of f{z) equals the type of exp{az'^)a{z). In the 
both cases the type of f{z) is greater or equal □ 

Remark C.4. If A is a quadratic or hexagonal lattice then z/ = and, consequently, 
a{z) = cr{z). Actually, in the former case we can suppose that Ui > 0, uj2 = iuJi- Then 
r]i = tt/uji, ri2 = —ixi/uJi [104, 18.14.8 and 18.14.10], hence i/ = 0. In the latter case 
we can suppose that uji = ke~'^'^^^, 102 = ke^'^^^, with A; > 0. Then 

27re*^/3 ^ 2Tre-'''/^ 

V^{uji + 002) ' V3(u;i + 002) 

18.13.16 and 18.13.19]. In this case, too, z/ = 0. 

Remark C.5. There exist lattices for which u and, consequently, a{z) 7^ cr{z). 
It suffices to consider a lattice with 772 = (such a lattice exists, see |104[ 18.3.10]). 
Then, by the Lagrange formula, 1771(^21 = 27r and hence |z/| = tt /2S. This means that 
the type of the a-function for such a lattice equals n/S. Since v depends on {(jJi,uj2) 
continuously any value of \v\ lying between and Ti/2S is realized by a convenient 
lattice. 
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